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ANALYTICAL  STUDIES  PERTAINING  TO  THE  MEASUREMENT  OF  LIGHT 


ABSORPTION  IN  AIR 

Bela  A.  Lengyel  and  R.  J.  Romagnoli 
San  Fernando  Valley  State  College 


ABSTRACT 


Experimental  methods  proposed  for  the  measurement  of  the  absorption 
coefficient  of  light  in  air  are  analyzed.  The  calculations  are  primarily 
intended  for  the  wavelength  regions  where,  for  certain  wavelengths,  the 
absorption  coefficient  to  be  measured  is  estimated  to  be  of  the  order 
of  10"9  cm-1.  The  effects  of  scattering,  heat  conduction,  static  and 
acoustic  pressure  variations  are  examined  and  are  related  to  the  problem 
of  detecting  the  signals  produced  by  the  absorption  of  heat  from  the 
laser  beam.  Both  pulsed  and  cw  schemes  are  considered.  The  measurement 
cf  the  static  pressure  rise  in  a  narrow  tube  appears  as  a  feasible,  al¬ 
though  marginal  methoo  which  offers  some  hope  of  success.  An  interfer¬ 
ometer  scheme  also  appears  feasible.  Finally,  the  linear  absorption 
coefficient  for  the  10. 6y  wavelength  COj  laser  beam  has  been  measured 


in  air. 
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1.  INTRODUCTION 


Propagation  of  light  at  very  high  intensities  produces  changes  in  the 
properties  of  the  medium  through  which  the  light  propagates.  These  changes 
in  the  medium  in  turn  affect  the  propagation  of  th«  beam.  Such  effects  have 
been  observed  in  liquids  and  gases  and  also  take  place  in  the  atmosphere  if 
laser  beams  of  sufficiently  high  intensity  were  caused  to  pass.3  The  laser 
beam  would  be  defocused  by  the  refractive  index  gradient  produced  in  the 
atmosphere  by  the  absorption  of  energy  from  the  beam.  The  work  of  Brueckner 
goes  beyond  this  first  order  effect  and  treats  the  combined  thermodynamic, 
mechanical,  and  optical  instability  resulting  from  the  interdependence  of 
density  distribution,  temperature  and  ray  path.  In  this  report  concern  is 
primarily  with  the  linear  tbsorption  coefficient  of  the  atmosphere  for  high 
intensity  laser  beams  of  various  frequencies,  and,  in  this  section,  with  the 
first  order  effect  of  defocusing.  Hence  a  review  of  the  principles  involved 
in  determining  beam  spread  is  presented  here. 

When  thermal  energy  is  deposited  in  the  atmosphere  (assumed  to  be  an 
ideal  gas)  by  absorption  from  a  laser  beam,  the  resulting  effects  on  pressure, 
temperature,  density,  and  molecular  velocity  are  determined  by:  (1)  the  iueal 
gas  laws,  (2)  the  equation  of  continuity,  (3)  Newton's  second  law  of  motion, 
and  (4)  conservation  of  energy.  These  effects  are  derived  in  many  stardard 
texts. k  applying  the  law  of  conservation  of  energy,  one  obtains  the  energy 

a.  K.A.  Brueckner:  I.D.A.  Laser  Summer  Study,  1963,  or  I.b.A.  Research 
Paper  P-42.  Report  on  Laser  Summer  Study,  The  Inst,  for  Radiation  Physics 
and  Aerodynamics.  U.C.S.D.,  Aug.  1965. 

b.  Knudsen  6  Katz:  FLUID  DYNAMICS  AND  HEAT  TRANSFER,  McGraw  Hill,  1958. 
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balance  equation  for  ideal  gases  (Eq.  2-55  of  Knudsen  and  Katz):*3 


PC'  £L  =  K72T  ♦  ♦  q'  ♦  (1.1) 

P  dt  dt 


where  p  =  mass  uensity 

T  =  absolute  temperature 

=  heat  capacity  at  constant  volume 
Cp  =  heat  capacity  at  constant  pressure 
t  =  time 
P  =  pressure 
Y  =  C£/C; 

K  =  thermal  conductivity 

4>  -  viscous  dissipation  function 

q'  =  rate  of  heat  generation  in  gas  per  unit  volume. 

Assuming  K  and  *  negligible,  and 
q'  r  ol, 

where  .  o  =  linear  absorption  coefficient  of  the  atmosphere  (ideal  gas) 
for  the  laser  beam, 

I  =  intensity  (power  per  unit  area)  of  the  laser  beam, 
then  equation  (1,1)  becomes. 


pC '  —  -  —  =  q ' 
p  P  dt  dt  M 


=  ol. 


For  an  ideal  gas 


C ' 
LP 


c  • 
<-v 


P 

PT 


(1.2) 


(1.3) 


Hence, 

cv(Y -!)<£  =  !<£  P  dp  . 

v  dt  p  dt  p2  dt 


(1.4) 
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Substituting  (1.4)  into  (1.2),  one  obtains  a  relationship  independent  of 
temperature  as  follows: 


, .  ( _ 1 _ [l  dP_  _  P  dpi  ]  _  d_P  _ 

P((C;(Y-1)  [p  dt  *  7  dtj  j  dt 


_i_  (V  dP  .  (,,_!)  dP  _  yP  do_)  =  aI 

Y-l  Y  dt  '  dt  p  dt ' 


(1.5) 


dP.lIdp  =  (y-l)  al. 
dt  p  dt 


Now  make  use  of  a  relationship  for  the  speed  of  sound  in  a  gas. 


C  j  -  YP/P, 


(1.6) 


and  one  obtains: 


(P  -  c\p)  =  (y-1)  al. 
dt 


(1.7) 


A  concise  treatment  of  the  effects  of  density  change  on  beam  spread 
and  instability  has  been  developed  by  Brueckner.a  His  basic  equations 
are  obtainable  from  (1.7)  by  taking  the  Laplacian  of  both  sides: 


d  (v2p  .  c2V2p)  =  (y-1)  V2Ia. 
at  ° 


(1.8) 


Now,  using  the  equation  of  motion 

-vp  .P5, 

dt 


(1.9) 


the  equation  of  continuity. 


-  H  •  V-(ov). 


(1.10) 


along  with  the  operator  identity, 


d_  =  i_  +  vV, 

dt  3t - 


(1.11) 
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where  v  is  the  fluid  velocity,  one  obtains  an  expression  independent  of 
pressure  as  follows : 


72p  =  -  v.(o^)  =  -  +  7 •  (v  &-) 

dt  dt  3t  dt 

=  t  v*V  —  Jt  7*[v  (—  +  v-7p)1 

at7  -  -  at  -  -  at  -  - 

=  — -  p(7-v)2  -  v7[2p7*v  +  v>7p]. 
oX.  “ 


(1.12) 


For  small  pressure  and  density  fluctuations,  neglecting  gradients  in  v, 
one  can  linearize  equation  (1.12)  by  the  approximation. 


72P 


Hence  the  equation  for  density  change  in  a  gas  traversed  by  a  laser  beam 
is: 

d?  (3t2  "  co*2>P  s  ^ Y”1 )  ^2Ia.  (2.13) 


If  density  changes  occur  slowly,  so  that 

32p 


<<  c272p , 


they  may  be  approximated  to  occur  at  constant  pressure,  and  equations  (i.7) 
or  (1.13)  simplify  to  the  following: 


=  -  2^.  la.  (1.14) 

3t 

Solutions  of  (1.13)  and  (1.14)  combined  with  the  eikonal  equation  have  been 
obtained  bv  Brueckner.3  Suffice  it  here  to  obtain  a  very  elementary  solu¬ 
tion  to  (1.14)  to  show  how  an  initial  first  order  effect  of  the  spreading 
of  an  initially  parallel  laser  beam  depends  upon  a. 
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Before  introducing  geometrical  optics,  it  may  be  helpful  to  arrive  at 
equation  (1.14)  from  even  more  elementary  considerations.  At  constant  pres¬ 
sure  the  heat  per  unit  volume  absorbed  from  the  laser  beam  by  an  ideal  gas  is 


32  !  al  dt  s  -CP  dT  =  pCP  dT‘ 

V  V  M 

where  N  =  number  of  moles 
M  =  mass  per  mole 

Cp  =  molar  specific  heat  at  constant  pressure 
Cv  =  molar  specific  heat  ot  constant  volume 
V  =  volume . 

Therefore, 

cl  *  ?E  p  i  T  ££  =  -  fL  dP 

M  M  dt  M  dt  y-1  dt 


(1.15) 


(1.16) 


where  use  is  made  of  the  expression  found  in  elementary  texts  for  the  speed 
cf  sound  in  a  gas: 


c 


2 

o 


=  t  =  Se  (c_ 
M  Cy  P 


V  s  ■  CP(T'1) 


T 

M’ 


and,  at  constant  pressure. 


dT  _  dp 
T  ‘  p" 


Equation  (1.16)  is  another  form  of  (1.14): 


3p  . 
3t  " 


al. 


(1.14) 


The  angular  ray  deflection  may  be  approximated  from  Snell's  law,  where 
L  is  the  coordinate  in  the  direction  of  propagation,  and  r  is  a  tx'ansverse 
coordinate,  as  illustrated  in  Figure  1.1. 
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Differentiating  Snell's  law  with  respect  to  L,  one  obtains: 


0  =  -  n  sin  0  1®.  cos  0 
3L  3L 


or, 


39  .  1  3n  a  .  i  3n  dL  .  1  3n 
3L  n  3L  n  3L  dr  n  ' 


(1.18) 


Variations  in  refractive  index  of  air  may  be  closely  approximated  by  equation 
(1.19):c 


fin  =  —  (n  -1). 
°o  0 


(1.19) 


Hence  (1.18)  becomes,  since  n  *  1* 


36  1  nQ-l  3 0  no-l  3p 


6L  n  o  3r 
o 


°o  3r 


(1.20) 


Now  equation  (1.14)  may  be  integrated  and  substituted  into  (1.20): 

0  -  PQ  =  -  alt, 

co 

where  I  is  the  time  average  intensity. 

Hence, 


30  _  (no_1)  ,  nat  31 


3L 


c2  3r 
co 


(1.21) 


The  change  in  intensity  due  to  spreading  of  the  beam  may  be  approximated  by: 


3I_  _  3__  j 


ro2 


3r  3r  0  [r(L)]2 
Hence  equation  (1.21)  becomes 


2I0  ro2 

r"r(L)Y3 


(1.22) 


36  _  2(n0-l)(Y-l)aI0t  r2  x 

3L  Pc  Cq 


[r(L))3  fr(L)] 


(1.23) 


c.  (Kuiper  and  Middlehurst:  TELESCOPES,  U.  Chicago  Press.) 


* 

l 


where  B  is  the  constant  so  defined: 

2(n0-l)(  Y-l)al"0t  r2 


B 


P0  <=o 


(1.24) 


Lquation  (1.23)  may  be  integrated  easily,  for  n  beam  initially  parallel, 
to  determine  the  beam  spread.  Prom  Fig.  1.1, 

tan  6  =  ^  s  6,  or  dr  *  6  dL. 
dL 

Substituting  these  relationships  into  11.23),  one  obtains: 


B 


_  B  dr 


de  =  —  dL  =  —  —  . 
r3  r3  6 

Then,  upon  integration,  one  obtains: 


Hen-e, 


1  e2  =  /*  ede  =  b  f  ^  .  I  ri-  -  i_- 

2  O  -o  *  i-0 


=  £  r—  — i 

ro  r3  2  ri  r2 


dr  _  „  _  /r2-ri  , 

dL  6 - 77T~ 


which  may  be  integrated  to  give: 

r 


I  =  5*  J  dL  . 

'9  'r  /  0  9  _  n 


L 


ro  /r2-r2  r0  o 


(1.25) 


Hence  it  follows  that  the  beam  radius  increases  with  distance  according  to 
to: 


r  =  r0  [1 


,  Jili* 


(1.26) 


and  the  beam  deflection  or  angle  of  spread  is: 


e  L 

e  =  /  dfl  =  /  B 


r3  L 

dL  =  -2.  / 

o  or3  ° 


dL 


.  B1* 


or 


r  r41fc  rn  T  r4!1* 

L2+  —2.  **  °  L2+  ^ 

B  L  B  . 


e  =  5k  [i  + 


(1.27) 
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Hence  it  is  seen  that  to  a  first  approximation  the  beam  spread  is  propor¬ 
tional  to  B^,  and  B  is  proportional  to  beam  power,  time,  and  a,  the  linear 
absorption  coefficient.  It  is  therefore  important  that  methods  be  avail¬ 
able  to  measure  values  of  a  for  various  laser  frequencies  where  high  in¬ 
tensities  are  physically  realizable. 


2a.  Statement  of  the  Purpose  of  the  Report. 


It  is  of  interest  to  know  the  amount  of  heat  evolved  in  a  material 
through  which  a  laser  beam  is  passing.  A  parallel  beam  of  light  is  attenu¬ 
ated  as  it  passes  through  matter  and  this  attenuation  is  mainly  attributed 
to  two  effects:  scattering  and  absorption.  Scattering  on  material  parti¬ 
cles  results  in  the  directional  redistribution  of  the  energy  of  the  incident 
electromagnetic  wave.  Absorption  results  in  the  conversion  of  electromag¬ 
netic  energy  into  other  forms  of  energy.  The  energy  lost  by  the  electro¬ 
magnetic  field  is  ultimately  converted  into  heat.  This  means  that,  possibly 
through  the  intervention  of  processes  of  atomic  and  molecular  excitation, 
the  energy  abstracted  from  the  electromagnetic  field  is  redistributed  among 
the  available  degrees  of  freedom  of  the  material  (gas),  so  that  it  is  ulti¬ 
mately  ur.distinguishable  from  energy  obtainable  from  contact  with  a  heat 
reservoir.  We  shall  assume  that  the  atomic  and  molecular  processes  res¬ 
ponsible  for  the  redistribution  of  energy  among  the  available  degrees  of 
freedom  take  place  so  fast  that  at  the  scale  at  which  we  contemplate  our 
observations  the  conversion  of  the  absorbed  energy  into  heat  may  be  regarded 
as  instantaneous.  It  should  be  noted  that  we  are  talking  here  about  energy 
redistribution  locally  and  not  about  its  redistribution  in  space.  The  latter 
generally  proceeds  at  a  much  slower  rate.  We  shall  thus  proceed  under  the 
assumption  that  the  energy  absorbed  from  the  electromagnetic  field  instantly 
appears  at  the  site  of  absorption  in  the  form  of  heat.  Such  an  assumption  is 
probably  valid  if  our  time  resoltuion  is  no  finer  than  10“7  sec. 

The  absorption  of  light  in  solids  and  liquids  is  usually  significant 
enough  so  that  a  direct  measurement  of  the  absorption  is  not  difficult.  In 
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gases,  such  as  pure  air,  the  measurement  of  the  absorption  of  light  is  very 
difficult,  because  the  rate  of  absorption  is  small  and  in  the  visible  region 
it  is  completely  masked  by  the  much  larger  effect  of  scattering.  The  situ¬ 
ation  is  somewhat  more  favorable  in  the  infrared  beyond  1.2  microns. 

The  optical  effects  of  scattering  and  absorption  combine  so  that  the 
intensity  of  a  plane  wave  passing  through  matter  is  described  by  the  equation 

Hx>  --  ,  (2.1) 

where  the  attenuation  coefficient  S  is  the  sum  of  two  terms: 

fi  *  a  +  a.  (2.2) 

Here  a  is  the  absorption  coefficient,  o  the  scattering  coefficient  and  x  is 
the  direction  of  propagation. 

In  the  spectral  region  where  the  ratio  a/a  is  large,  it  is  not  feasible 
to  determine  a  by  measuring  the  attenuation  of  a  light  bea>r.  through  matter. 
The  value  of  a  must  thus  be  determined  not  by  the  change  produced  in  the  ra¬ 
diation  but  by  the  change  produced  in  the  absorbing  matter.  We  already  noted 
that  the  absorption  of  radiation  results  ivi  an  increase  in  the  thermal  energy 
of  the  absorbing  matter,  therefore  we  shall  address  ourselves  to  problems 
pertinent  to  the  detection  and  measurement  of  small  quantities  of  thermal 
energy  deposited  in  the  path  of  a  narrow  light  beam  of  the  type  that  may  be 
obtained  from  a  laser. 

It  has  been  estimated  that  for  certain  wavelengths  in  and  around  the 
visible  region  of  the  spectrum  a  is  of  the  order  of  10“9  cm-1.  Scattering 
in  pure  air,  as  calculated  from  Rayleigh's  theory,  produces  a  much  higher 
attenuation.  For  green  light  (0.5  microns),  o  =  1.7  x  10" //cm  for  air  free 
of  all  contaminants.  We  reserve  the  quantitative  data  for  the  next  section, 
but  we  note  here,  that  around  10  microns  the  situation  is  quit*  different. 
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First,  significant  and  spectrally  selective  absorption  is  likely  to  be  en¬ 
countered  because  of  the  presence  of  the  molecular  bands  of  the  constituents 
of  air  in  this  region.  Second,  Rayleigh  scattering  on  molecules,  which 
varies  with  the  fourth  power  of  the  wavelength,  is  reduced  to  negligible 
values  around  10  microns. 

This  report  contains  calculations  pertaining  to  thermal  and  acoustic 
effects  expected  as  a  result  of  the  passage  of  a  powerful  laser  beam  throupji 
air.  The  configurations  examined  are  such  which  are  likely  to  be  chosen  in 
experiments  aimed  at  the  determination  of  i  in  the  0.5  to  1.2  micron  spectral 
region.  We  have  included  here  the  analysis  of  the  pressure  change  in  a  com¬ 
pletely  enclosed  cylinder  whose  walls  are  kept  at  a  constant  temperature,  the 
calculation  of  the  rate  of  dissipation  of  heat  around  a  laser  beam  by  means 
of  thermal  conduction,  the  generation  of  a  quasi-static  pressure  increment, 
and  an  acoustic  pulse  by  a  laser  pulse  and  some  problems  pertaining  to  the 
detection  of  an  acoustic  pulse  or  a  sequence  of  acoustic  pulses.  The  material 
as  assembled  here  does  not  represent  a  complete  investigation.  It  is  a  group 
of  exploratory  calculations  preliminary  to  the  experimental  determination  of 
the  absorption  coefficient  of  light  in  air. 

2b.  Relevant  Properties  of  Air. 

In  this  section  we  have  collected  such  physical  data  concerning  air  that 
seem  relevant  for  computations  in  later  sections,  or  that  may  be  needed  in 
evaluating  an  experimental  scheme  proposed  for  the  measurement  of  the  absorption 
coefficient  a. 

Atmospheric  air  free  of  dust  and  other  contaminations  is  a  '.dxture  of 
gases,  primarily  N2,  02,  Ar,  C02  and  water  vapor.  Other  gases  are  present  in 
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such  smII  proportions  that  thair  praaanca  dots  not  noticeably  affect  the 
physical  properties  of  air  that  are  of  concern  to  us.  Of  the  gases  listed 
here  only  the  proportion  of  water  vapor  is  variable. 

In  what  immediately  follows  we  shall  assume  that  we  are  dealing  with 
clean,  dry  atmospheric  air.  By  this  substance  we  mean  a  gas  mixture  from 
which  all  water  vapor  has  been  removed  and  which  contains  the  gases  N2,  02, 

A r  and  C02  in  the  proportion  as  they  are  universally  found  in  the  atmosphere1. 
We  note  in  particular  that  such  air  contains  0.033  percent  C02  by  volume. 
Although  this  proportion  seems  small,  it  is  well  to  remember  that  a  relatively 
small  amount  of  CO^  or  water  vapor  significantly  affects  the  optical  properties 
of  air. 

It  is  customary  to  specify  the  physical  properties  of  gases  at  standard 
temperature  and  pressure;  i.e.  at  0°C  and  760  torr.  At  these  conditions  the 
properties  of  dry  atmospheric  air  relevant  to  this  study  are  as  follows: 
Density:  pQ  *  1.293  x  10"3g  cm"3. 

Specific  heats:  Cp  *  0.240,  cv  *  0.172  cal  g"1  °C"1, 

Sound  velocity:  cQ  *  3.31  x  104  cm  sec"1. 

Heat  conductivity:  K  *  5.61  x  10”5  cal  cm"1  sec"1  °C"1, 

Index  of  refraction  for  A  ■  5000  X:  n  *  1.000294. 

At  standard  temperature  and  pressure  the  particle  density  of  every  ideal 
gas  is  L  *  2.69  x  10 13  cm-3.  (Loschmid's  number.)  The  variation  of  n  -  1 
with  wavelength  (dispersion)  has  been  studied  and  is  tabulated  for  a  tempera¬ 
ture  of  15°C,  which  must  be  a  convenient  laboratory  temperature  in  Sweden. 

The  following  data  applicable  at  15°C  are  obtained  from  the  work  of  Edlln:2 

T.  Handbook  of  Chemistry  and  Physics,  46th  ed.  p.  F-116. 

2.  B.  Edlln,  The  dispersion  of  standard  air,  J.  Opt.  Soc.  Amer.  43_,  339,  1953. 
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X(Angstroms)  5000  6000  7000  8000  9000  10000 

(n  -  1)  x  104  2.78  2.76  2.75  2.75  2.74  2.74 

At  the  more  comfortable  laboratory  temperature  of  20°C  the  following 
data  are  applicable  for  760  torr  pressure: 

Density:  p2q  =  1.205  g  cm-3, 

Sound  velocity:  c2Q  =  3.44  x  104  cm  sec-1, 

Loschmid's  number  L  =  2.51  x  i0lg  cm-3, 

The  value  of  n  -  1  varies  from  2.73  x  10“4  at  5000  A  to  2.69  x  10“  at  10000 

The  attenuation  of  light  caused  by  scattering  on  the  molecules  of  air  is 
calculable  for  Rayleigh's  formula3 

o  =  .  (2.3) 

3LX4 

At  standard  temperature  and  pressure  we  obtain  the  following  values  of  o: 
X( Angstroms)  5000  6000  7000  8000  9000  10000 

o  x  107(cm" 1 )  1.70  0.81  0.44  0.25  0.16  0.10 

Both  (n-1)  and  L  are  proportional  to  p,  therefore  o  is  also  proportional 
to  p .  At  constant  pressure  o  is  inversely  proportional  to  the  absolute  tem¬ 
perature.  Therefore  the  values  of  o  at  t  =  20°C  and  760  torr  are  about  7 
percent  smaller  than  the  above  tabulated  values. 

3.  Absorption  Experiment  With  a  Steady  Laser  Beam. 

We  shall  analyze  the  problems  involved  in  the  measurement  of  a  by  means 
of  the  following  experiment:  A  certain  volume  of  air  is  confined  at  a  known 

T.  J.W.  Strutt,  Lord  Rayleigh,  On  the  transmission  of  light  through  an 
atmosphere,  Phil.  Mag.  47^  374,  1899;  Scient.  Papers  4.,  397. 
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initial  pressure.  A  laser  beam  passes  through  the  confined  volume  and  heats 
the  enclosed  air  at  a  rate  proportional  to  a.  In  a  steady  state,  the  heat 
losses  of  the  system  balance  the  heat  gained  from  the  laser.  The  increase 
of  the  static  pressure  within  the  vessel  is  measured  and  is  used  for  computing 
the  rate  at  which  heat  is  added  to  the  air. 

The  simplest  experimental  embodiment  of  this  idea  consists  of  a  cylin¬ 
drical  tube  kept  at  a  constant  temperature.  The  laser  beam  passes  along 
the  axis  of  the  cylinder  heating  a  smaller  cylindrical  volume  of  air  con¬ 
centric  with  the  confining  tube.  A  steady  state  is  established  when  the 
rate  at  which  heat  is  conducted  to  the  walls  is  equal  to  the  rate  at  which 
heat  is  absorbed  from  the  laser  beam. 

In  order  that  a  technique  of  this  type  yield  a  valid  measurement  of 
a,  certain  physical  conditions  must  be  satisfied  which  are  not  easy  to 
achieve  experimentally.  First,  the  confining  tube  must  be  long  enough  so 
that  special  thermal  effects  at  the  ends  of  the  tube  do  not  significantly 
affect  the  pressure  of  the  confined  air.  Laser  light  is  brought  in  through 
solid  windows  whose  density  is  about  1000  times  that  of  the  confined  gas. 

It  is  not  unreasonable  to  assume  that  absorption  per  unit  volume  in  the 
windows  will  be  about  1000  times  as  great  as  absorption  in  the  gas. 

Second,  for  A  <  1.2y  more  laser  light  will  be  scattered  than  directly  ab¬ 
sorbed.  The  energy  of  the  scattered  radiation  must  be  removed,  it  must 
not  be  allowed  to  add  to  the  thermal  energy  of  the  confined  gas.  This 
removal  could  be  accomplished,  for  example,  by  making  the  inner  surface 
of  the  tube  black  and  by  making  the  tube  of  a  material  with  good  thermal 
conductivity. 
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Assuming  that  these  experimental  difficulties  are  mastered,  we  calcu¬ 
late  the  pressure  increase  dp  attained  in  the  vessel  in  the  steady  state 
as  a  result  of  energy  being  absorbed  from  the  laser  beam  at  a  constant 
rate  and  conducted  away  through  the  air  to  the  walls  of  the  vessel  which 
are  kept  at  a  constant  temperature  T  .  We  assume  a  uniform  circular  laser 
beam  of  radius  a  in  a  vessel  of  radius  b.  Let  the  power  incident  in  the 
laser  beam  be  W,  then  the  intensity  of  the  beam  is  I  =  W/a2*.  Heat  is 
deposited  in  the  unit  volume  of  the  gas  at  the  rate  aJI,  where  J  =  0.230 
cal/ joule. 


Once  the  steady  state  is  reached,  no  mechanical  work  is  done.  There¬ 
fore  the  heat  balance  equation  is  obtained  by  equating  the  rate  at  which 
energy  is  absorbed  to  the  rate  at  which  it  is  conducted  away.  (It  is 
assumed  that  convection  is  neglected.)  The  rate  of  heat  generation  per 
unit  volume  is 

q'(r)  =  aJI  when  r  <  a, 

=  0  when  r  >  a . 

The  heat  balance  equation  for  a  volume  V  is  then 

Jq'dv  =  -  K  ^VT-ndS,  (3.1) 

V  S 


where  dS  is  the  surface  element,  n  is  the  unit  normal  to  the  surface,  and 
K  the  thermal  conductivity.  It  follows  then  in  a  well-known  manner  that 
V2T  =  -  q'/K.  (3.2) 

Since  the  only  relevant  coordinate  in  this  problem  is  r,  the  distance 
from  the  cylinder  axis,  we  have 


1  d_  dT_  _  q 1  ( r ) 
r  dr  r  dr  K 


C  3 . 3  ) 
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This  differential  equation  is  to  be  solved  with  the  boundary  condition 
T  =  Tq  for  r  -  b  and  the  solution  is  required  to  remain  finite  for  r  <  b. 
Although  the  solution  is  easily  found,  we  do  not  need  its  explicit  form 
because  we  are  only  interested  in  the  pressure  rise  in  the  vessel  and  this 
is  proportional  to  the  increase  of  the  average  temperature. 


In  place  of  T  we  introduce  the  dimensionless  variable 

v(r)  s  sr(T '  T=>- 


(3.4) 


(3.5) 


Equation  (3.3)  is  th"n  replaced  by 

Id  dv  .  . 

r  ^r  dr 

where  g(r)  =  1  for  r  <  a  and  0  for  r  >  a.  The  boundary  condition  is 
v(b)  =  0.  Then  we  obtain  by  integration  of  (3.5) 

r 


r  ~  =  -  /  gfr^r'dr*  =  -  f(r), 
dr  o 

where  f(r)  =  ^  r2  for  r  <  a,  and  f(r)  =  a2  for  r  >  a. 


(3.6) 


The  average  increase  in  gas  temperature  over  the  entire  vessel  is 
proportional  to 

b  2ir  L  b 

v  =  -  j  If  v(r)rdrd6dz  =  — j-  J  v(r)rdr.  (3.7) 

o  o  o  b  o 

Integrating  by  parts  ard  noting  that  v(b)  =  0,  we  get 

b  b  b 

I  v(r)rdr  =  -r.  f  r2  IX.  dr  =  i  /  f(r)rdr. 
o  2  o  dr  2  o 


U  a  u  2  2  2 

Since  /  f(r)rdr  =  fir3  dr  ♦  /  ia  rdr  =  a  b  [2  -  a  ]t  we  have 

n  O  2  '2  8b7 


.2  -2 

v  =  Cl  -  5—]. 

4  2b2 


(3.8) 


Therefore  the  average  increase  o  '  temperature  in  the  vessel  is 

dT  =  [1  -  — ].  (3.9) 

4K  2b2 

The  final  expression  for  the  increase  of  pressure  in  the  chamber  is  ob¬ 
tained  from  (3.9)  by  using  the  gas  law  and  by  reintroducing  the  incident 
power  W  =  a2ul. 

dp  = - (1 - *-).  (3.10) 

Tc  4itK  2b2 

Ordinarily  the  expression  in  the  parentneses  is  nearly  equal  to  one. 
Assuming  a  laser  with  100  watts  useful  output,  an  initial  temperature 
T^  =  300°K  and  a  =  10-9,  the  expression  in  front  of  the  parentheses  is 
1.1  x  10- 7  atmospheres .  A  pressure  change  of  this  order  of  magnitude 
should  be  detectable  provided  one  can  keep  the  temperature  of  the  walls 
constant  to  within  0.3  x  10-1*  °C .  It  might  be  quite  difficult  to 
accomplish  this  for  a  tube  as  long  as  1  meter. 

4.  Dissipation  of  Heat  in  Air  by  Conduction. 

In  preparation  for  the  discussion  of  experiments  involving  the  heat¬ 
ing  of  air  columns  by  means  of  short  laser  pulses  we  shall  calculate  the 
rate  at  which  thermal  energy  initially  communicated  to  an  air  column  is 
dissipated  by  thermal  conduction  alone. 

The  calculation  is  simplified  by  neglecting  for  the  purpose  of  this 
calculation  the  aerodynamic  properties  of  air.  Thus  we  shall  be  concerned 
with  a  solid  substance  characterized  by  a  heat  conductivity  K,  density  p. 
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and  specific  heat  cv. 


It  is  shown  in  textbooks4  that  the  temperature  v(x,y,z,t)  satisfies 
the  differential  equation 


72v  -  I  * 
K  3t 


0, 


(4.1) 


provided  that  there  are  no  sources  or  sinks  of  heat  in  the  region  under 
consideration.  Here  k  =  K/pcv  is  the  diffusivity  of  the  substance. 


We  will  assume  a  geometry  of  cylindrical  symmetry  and  calculate  the 
cooling  of  ar.  infinite  cylinder  of  radius  a,  whose  temperature  at  the  time 
t  =  0  is  one  degree  above  the  temperature  of  its  surroundings.  All  space 
is  assumed  to  consist  of  the  same  material  characterized  by  the  constant 
k.  The  problem  then  requires  the  use  of  only  one  spatial  variable,  r, 
the  cylindrical  polar  coordinate.  The  differential  equation  becomes 


and 


the 


initial  condition 
f (r)  =  1 
=  0 


3r2  r  3r  / 
is  v(r,0)  =  f(r), 
for  r  <  a, 
for  r  >  a . 


where 


(4.2) 


As  lonp  3  no  walls  are  present,  no  boundary  conditions  are  required. 
This  differential  equation  is  readily  solved  by  the  Hankel  transformation 
technique5  and  the  following  result  is  obtained: 

*T!  H.S.  Carslaw  and  J.C.  Jaeger,  Conduction  of  Heat  in  Solids,  Oxford,  1569, 


2nd  ed.  p. 

9. 

Ref.  4,  p. 

460. 
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(4.3) 


v(r,t)  =  /  F0(o  )oJ0(or)e_lCO  tdo, 
o 


where 


r0(o)  =  /  rJQ(or)f ( r)dr . 
o 


(4.4) 


In  the  present  case 


F0(o)  =  /  rJ0(or)dr  =  —  Jj^(oa). 
o  a 


(4.5) 


Hence  the  expression  for  the  temperature  distribution  is 


r(r,t)  =  a  /  J^oa  )JQ(  or)e  "’co  x  do. 


(4.6) 


It  is  interesting  to  calculate  the  variation  of  temperature  at  the  axis. 
From  (4.6)  we  have  for  r  =  0 


v( 0 ,  t )  =  a  /  J0(oa)e-ico  1  do. 
o 

This  expression  is  changed  by  integration  by  parts  into 


(4.7) 


r(0,t)  =  1  -  2<t  /  J0(oa)oc"'co  x  do. 
o 


From  the  standard  formula6 


(4.8) 


/  J  (oa)e"P  0  odo  =  -i—  e  UP 
o  °  2p2 


(4.4) 


we  ge  t 


v(0,t)  =  1  -  e 


4ic  t 


(4.10) 


6.  G.N.  Watson,  Theory  of  Bessel  Functions,  Cambridge,  University  Press, 
1958,  p.  393. 
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Thus  the  temperature  at  the  center  line  decreases  with  t  according  to  tnc 
formula  1  -  e-1/*,  where  t  =  a2/4ic.  The  quantity  7  has  the  dimension  of 
time.  It  is  the  time  during  which  the  temperature  decreases  by  the 
amount  1/e.  Although  we  are  not  dealing  here  with  a  simple  exponenti.il 
decay  of  temperature,  it  is  still  reasonable  to  regard  t  as  a  thermal  r> 
taxation  time.  It  is  proportional  to  the  cross  sectional  area  of  the 
beam . 


numerical  computation  of  the  thermal  relaxation  time  requires  the 
knowlo Ige  of  <  for  air  in  the  applicable  temperature  ranee.  This  quan¬ 
tity  entered  first  in  equation  (4.1)  which  is  taken  from  the  theory  ot 
heat  conduction  applicable  to  solids.  When  this  theory  is  applied  to 
a  pas ,  the  question  arises  what  kind  of  specific  heat  is  to  be  used  in 
the  equation 

K  =  --  (4.11) 

pc 

flint  serves  as  a  definition  of  <2.  Since  we  are  contemplating  a  steady 
state  situation  in  which  no  mechanical  motion  taker,  place,  no  work  is 
being  done  and  the  entire  heat  energy  that  enters  a  volume  of  pas  is 
converted  into  internal  energy.  The  specific  heat  at  constant  volume 
is  the  one  which  is  applicable  to  such  a  situation.  For  air  at  standard 
pressure  and  temperature,  <  =  0.255  cm2  sec-*.  Although  depends  on  the 
tomoerature ,  this  value  of  k  is  probably  a  fair  approximation  to  the 
correct  value  applicable  at  20°C.  With  the  value  tc  =  0.255  cm2  sec  1 
we  obtain  the  following  approximate  values  of  relaxation  times  and  associ 
ated  frequencies  applicable  to  beams  of  convenient  cross  sections. 
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beam  diarr.  (2a)  mm  0.4 

0.6 

0.8 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 

<elax.  time  (t)  msec  0.392 

0.882 

1.57 

2.45 

3.53 

4 . 80 

6.27 

7.94 

l.rt 

frequency  (t_1)  Hr  .’SsO 

1134 

637 

408 

283 

208 

160 

1  ?r 

10. 

In  a  dynamic  situation. 

when 

heat  is 

added 

to  a 

column 

of  ga 

s ,  not 

only  heat  conduction  tak«s  place  hut  the  pas  is  set  in  motion.  This  trans¬ 
ient  phenomenon  was  not  taken  into  account  so  far.  In  what  follows  later, 
the  dynamic  •  i.enom^non  is  of  great  importance  and  the  results  of  this  section 
arc.'  used  for  guidance  concerning  the  situations  when  the  effect  of  heat  ^in¬ 
duction  may  be  neglected.  Clearly  one  can  onLy  justify  neglectinp  heat  con¬ 
duction  it  the  dynamic  phenomenon  (oscillation)  occurs  so  fast  that  only 
negligible  cooling  takes  place  during  one  cycle.  Cucn  will  be  the  case 
when  the  frequency  of  the  sound  wave  generated  exceeds  the  appropriate 
c.hai  acteristic  frequency  t"1  =  ^K/a7 . 

When  che  air  colurn  is  not  free,  as  heretofore  assumed,  but  confine:  is 
a  pipe  kept  at  a  constant  temperature?,  the  tie t hod  of  calculation  is  different. 
The  differential  equation  (4.1)  is  applicable  in  this  case,  but  a  boundary 
condition,  absent  in  the  case  of  the  free  air  column,  must  lie  added,  for  a 
cylindrical  vessel  of  inner  radius  b,  the  boundary  condition  is  v(b,t)  =  0, 
provided  that  the  vessel  is  kept  at  a  constant  temperature  T0 .  (Mere  v 
denotes  the  excess  over  this  temperature.)  The  solution  of  the  heat  con- 
duction  problem  for  a  domain  with  a  boundary  is  well  known  .  It  is  obtained 
by  expanding  the  function  to  be  determined  in  terms  of  the  eigenfunctions 
of  the  differential  equation  (4.1)  with  the  given  boundary  conditions. 

In  the  case  of  heat  conduction  within  a  cylinder  of  radius  b  the  result  is 

7~  ReT.~:  "pp-  1M8- 199. 
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(4.12) 


v(r,t)  -  l  AnJ0(anr)e  n  , 
n  =  l 

where  the  numbers  an  are  determined  from  the  boundary  condition  J.  (anb)  =  0, 
and  the  coeff ic j.ent j  Af)  are  calculated  from  the  exDansion 


f(r)  =  I  An,Io^nr),  (4.13) 

n=l 


which  the  explicit  form  of  the  initial  condition.  The  A  are  deter¬ 
mined  by  the  routine  procedure  applicable  to  expansions  accordinp  to  or¬ 
thogonal  functions.  Their  values  depend  on  a  and  b.  In  the  pcnernl  case 
the  value  of  An  does  not  drop  off  rapidly  with  increasinp  n.  We  must 
therefore  take  into  account  a  large  number  of  terms  in  (4.1?).  Each  term 
decays  in  time  accordinp  to  a  different  exponential  rate  The  t itne* 

constants  associated  with  these  modes  are  -^  =  (ca^)*1,  and  an  is  deter¬ 
mined  from  the  equation 

rtn;’  =  Xn,  ( *’  * lu  > 

where  x  is  the  n-th  root  of  the  equation  JQ(x)  =  0.  The  values  of  the 
first  five  x's  are:  2.405,  5.520,  8.654,  11.792  nd  14.ni3,8  The  time 
constants  are 


T 


n 


(4.15) 


We  note  that  the  role  of  the  tr's  is  not  the  same  as  that  of  t  in 
the  case  of  the  free  cylinder  because  in  the  present  case  the  time  de¬ 
pendence  of  a  mode  is  of  the  form  exp(-t/tn)  and  not  1  -  exn(-t/t). 

By  means  of  a  physical  argument  we  can  readily  infer  that  the  confined 

8.  A.  Gray  and  G.B.  Mathews,  Bessel  Functions,  2nd  ed.,  Dover  Publ.,  New 
York  1966,  p.  300. 
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cylinder  will  cool  more  rapidly  than  the  free  cylinder  and  that  the  re¬ 
laxation  time  calculated  for  the  free  cylinder  will  he  an  approximation 
of  the  effective  relaxation  time  of  a  confined  cylinder  for  b  >:■  a. 

5.  Pressure  Rise  in  a  Tube  Caused  bv  a  Single  Laser  Pui.se. 

One  night  atten,  t  to  confine  air  in  a  cylindrical  tube,  shoot  a  laser 
beam  through  and  measure  the  rise  in  pressure  in  the  ube.  The  expected 
increase  in  pressure  is  proportional  to  the  energy  deposited  in  the  air 
1) v  the  laser  beam.  Hence  one  expects  that  the  absorption  coefficient  a 
might  1<  ?  measured  in  such  an  experiment. 

It  must  be  remembered  that  the  walls  of  the  tube  holding  the  air  are 
solid  and  have  a  much  larger  heat  capacitance  than  the  confined  air.  So 
eventually  as  thermal  equilibrium  is  approacned ,  the  thermal  energy  de¬ 
posited  in  the  air  is  conducted  to  the  walls,  temperature  and  pressure  in 
the  vessel  return  to  their  values  at  the  beginning  of  the  experiment.  The 
approach  to  thermal  equilibrium  may  be  thought  of  as  consisting  of  a  fast 
process  cf  pressure  equalization  followed  by  a  slower  process  of  tempera¬ 
ture  equalization  by  heat  conduction.  The  pressure  equalization  m-’  be 
characterized  bv  the  time  constant  2b/c  where  b  is  the  radius  of  the  tube 
and  c  the  velocity  of  the  sound.  The  rate  at  which  heat  conduction  leads 
to  thermal  equilibrium  must  be  calculated  in  the  manner  described  in  Section 
4.  The  time  constant  t  =  L2/4k  certainly  reveals  the  correct  order  of  mag¬ 
nitude.  The  measurement  of  pressure  increase  should  take  place  after  the 
elapse  of  the  time  2b/c  following,  the  passage  of  the  laser  beam  but  consid¬ 
erably  before  the  elapse  of  the  time  t.  Tor  a  tube  of  1  cm  radius, 

2b/c  =  6  x  10" 5  sec  and  t  =  1  sec. 


The  laser  beam  will  generally  be  much  narrower  than  the  tube.  Let 
the  radius  of  the  beam  cross  section  be  a,  and  let  us  assume  that  we  are 
dealing  with  a  very  short  laser  pulse.  Then  at  the  end  of  the  pulse  we 
have  an  air  cylinder  of  radius  a  within  which  the  pressure  has  been  raised 
to  p  +  dp  and  the  temperature  to  T  +  dT,  while  outside  this  cylinder  the 
pressure  is  p  and  the  temperature  is  T.  According  to  the  standard  equation 
of  heat  balance 

cypoa2Trd7  =  aJQ,  (5.1) 

where  Q  is  the  energy  of  the  laser  pulse  in  joules,  J  is  the  conversion 
factor  0.239  cal/joule,  and  pQ  is  the  density  of  the  air.  From  the  ideal 
gas  law  and  equation  (5.1),  we  get 


aJQ 


Vo**** 


(5.2) 


On  substituting  the  appropriate  physica)  constants  of  air,  we  get 
for  a  =  0.1  cm  and  T  =  300°K: 

^  *  1.14  x  102  aQ.  (5.2) 

P 

Thus  a  laser  pulse  of  1  joule  and  an  a  =  10-9  cm-1  would  give  a  relative 
pressure  incrr >se  of  one  part  in  ten  million.  It  seems  then  that,  using 
a  powerful  laser  of  several  joules  output,  an  a  of  the  order  of  10-9  cm-1 
would  be  just  barely  detectable.  This  would  be  the  case  if  the  entire  con¬ 
fined  space  would  be  within  the  beam.  Actually  the  beam  diameter  must  be 
considerably  smaller  than  the  diameter  of  the  tube  that  keeps  the  air  con¬ 
fined.  Consequently  the  pressure  rise  observable  in  the  tube  is  much  less 
than  the  pressure  rise  in  the  cylinder  swei.t  out  by  the  beam. 
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An  estimate  of  the  pressure  rise  antit.ipated  in  the  tube  may  be  ob¬ 
tained  by  ascribing  the  entire  p,as  in  the  tube  some  sort  of  average  tem¬ 
perature  T.  The  heat  balance  equation  (5.1)  is  then  replaced  by 

cvP0')2l,dT  =  aJQ.  (5.4) 

This  equation  differs  from  (5.1)  insofar  as  b2  takes  now  the  place  of  a2. 
Consequently  the  pressure  rise  calculated  from  (5.4)  will  be  reduced  in 
proportion  to  the  ratio  of  the  cross  sectional  ratio  of  the  beam  to  the 
cross  sectional  area  of  the  tube.  If,  for  example,  a  =  0.1  cm  and  b  =  1.0 
the  valje  of  dp/p  is  reduced  by  a  factor  of  100,  thus  driving  it  completely 
beyond  reach  for  measurement . 


It  could  be  objected  that  the  above  estimate  of  dp/p  is  based  on  a 
dubious  assumption  of  an  average  temperature  and  the  appl ication  of  thermo¬ 
dynamics  to  non-equilibrium  states.  The  deduction  may  be  replaced  by  the 
following  one:  After  irradiation,  the  pressure  of  the  air  within  the  beam 
is  pj,  its  volume  is  v^.  In  the  rest  of  the  vessel  the  pressure  is  pQ . 

The  volume  of  the  rest  is  V2-  Both  volumes  change  during  the  process  of 
pressure  equalization  which  results  in  a  finai  pressure  P.  Assuming  no 
heat  exchanpe, 

vi  =  v!  and  p^//-'r  v_  =  P^Y  v'. 

11  1  ro  /  2 


Moreover,  the  total  volume  does  not  change, 
Hence 


)1/Y  .  Pi 


1/Y 


Po 


1/Y 


+  v- 


therefore  vj  +  vj 


V1  +  v2 • 


(5.5) 


and 


pl/Y 


V1 

V1  +  v2 


(5.6) 


cm 
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from  this  last  equation  it  also  follows  that  in  first  approximat ion 

P  -  p  ^  (p-p).  (5.7) 

°  vi  +  v2  1  ° 

6.  Generation  of  Sound  by  a  Laser  Pulse. 

A  laser  pulse  is  sent  through  a  pas  which  has  a  small  absorption  co¬ 
efficient  for  the  radiation  sent  through  it.  As  a  consequence  of  absorption 
a  column  of  gas  is  warmed  and  the  surrounding  gas  is  set  in  motion.  It  is 
desired  to  determine  the  sound  energy  generated  by  the  pulse  and  to  determine 
the  spectral  distribution  of  the  sound. 

The  solution  of  the  problem  is  expected  to  contain  the  physical  para¬ 
meters  of  the  gas  which  normally  determine  the  velocity  of  propagation  of 
sound  in  the  pas,  the  rate  at  which  radiated  energy  is  absorbed  in  the  pas, 
the  length  of  the  laser  pulse  and  the  cross-section  of  the  laser  beam. 

It  will  be  assumed  that  the  gas  is  an  ideal  pas,  that  the  laser  beam 
is  of  circular  cross-section,  that  the  rate  of  energy  dissipation  and  the 
pulse  length  are  such  that  the  resulting  temperature,  pressure  and  density 
changes  may  be  treated  as  differentials  of  the  first  order  so  that  the 
usual  approximations  employed  in  the  theory  of  sound  retain  their  validity. 

The  effects  of  heat  conduction  will  be  neglected:  a  simplification 
which  may  only  be  questionable  when  it  is  applied  to  low  frequencies.  The 
low  frequenc/  1  of  this  approximation  is  the  characteristic  frequency 

associated  with  the  coolinp  of  the  heated  column  by  conduction.  This  was 
calculated  in  Section  *« . 
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The  propagation  of  sound  waves  is  described  by  a  solution  of  the 
wave  equation 

-  c2?2$  =  0  (6.1) 

3t2 

which  is  satisfied  by  the  velocity  potential  $  as  well  as  by  the  pressure 
p  and  the  condensation  s  in  the  gas.  The  differential  equation  (6.1)  is  a 
homogeneous  one,  the  solution  #(x,y,z;t)  is  determined  if  one  knows  the 
value  of  <f>  and  9 <f> / 9 1  at  the  time  t  =  0. 


The  physical  problem  to  be  solved  here  must  be  governed  by  the  dif¬ 
ferential  equation  (6.1)  in  the  region  where  no  heating  takes  place.  In 
the  region  traversed  by  the  laser  beam  we  expect  to  obtain  an  inhomogeneous 
differential  equation  whose  solution  is  to  be  obtained  for  a  given  right 
side.  Moreover  the  solution  is  to  vanish  for  t  =  0;  it  should  represent 
an  outgoing  wave  away  from  the  region  traversed  by  the  laser  beam  and  is 
to  satisfy  the  usual  continuity  conditions. 


The  starting  point  in  acoustics  is  the  continuity  equation  for  the 

fluid 


+  V-pV  =  0 
3t 


(6.2) 


and  the  equation  of  motion  for  the  fluid  particles: 


P  =  -  ?p.  (6.3) 

at 

Here  the  vector  V  denotes  the  velocity  of  the  fluid,  p  its  pressure,  p 
its  density.  All  these  quantities  are  functions  of  position  and  time. 
It  is  customary  to  restrict  the  problem  to  small  velocities  and  small 
displacements  from  static  or  average  pressure  and  density  pQ  and  pQ  . 
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With  these  restrictions  equations  (6.2)  and  (6.3)  are  replaced  by 


and 


3p  — 

it  +  poV'V  =  0 


Po  ft-  ♦  Vp  -  0. 


we  introduce  the  velocity  potential  <J>  defined  by 


(  6 . 4 ) 


(6.5) 


V  = 


(6.6) 


then  (6.5)  becomes 


-P_  rr  +  VP 

O  ot 


=  0, 


(6.7) 


or 


v(|£  -  E_)  =  o. 

p0 


(6.8) 


We  write  p  =  pQ  +  p ' ,  where  pQ  is  a  constant  and  p'  is  a  variable  quantity 
small  compared  to  pQ.  Then 


"  P0 


(6.9) 


This  last  equation  relates  the  pressure  variation  to  the  velocity  potential. 

Now  the  continuity  equation  (6.4)  is  rewritten  by  introducing  the  con¬ 
densation  s  defined  as  follows 


p  =  pc ( 1  +  s). 


(6.10) 


or 


s  = 


P  -  Pc 
Po 


(6.11) 


Then  (6.4)  becomes 


ds  +  y.y 
3t 


=  0, 


(6.12) 


or 
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(6.13) 


3-  =  V26. 

at 


In  the  classical  case  the  equation  of  sound  propagation  ( F.q .  (6.1)) 
is  found  bv  assunir.r  that  the  pressure  variation  p'  is  proportional  to 
the  condensation  s.  The  factor  of  proportionality  is  found  by  thermo¬ 
dynamic  argument  to  be  equal  to  yp0.  With  this  result  one  gets 

p '  =  YP0s .  (6.14) 


lienee  from  (6.9) 


d2<t> 

\3 

to 

3s 

i 

Csi 

4— 1 

rc 

Po 

3t 

(6.16) 


finally 


a2$  _  Yp0 
3t2 


?2  i 


(6.16) 


which  is  identical  with  equation  (6.1).  Its  derivation  is  based  on  the 
assumption  that  no  thermal  exchange  takes  place. 


When  the  gas  is  Locally  heated  equation  (6.14)  is  no  longer  valid. 

We  must  calculate  the  change  in  pressure  when  both  the  volume  is  changed 
and  a  certain  amount  of  heat  is  added.  If  the  heat  is  aided  reversibly, 
we  write  dO  =  TdS.  Then 

dp  =  (l£)  dV  +  (|£-)  dS.  (6.17) 

3V  q  9S  y 

Equation  (6.14)  is  a  special  case  of  (6.17)  for  dS  =  0. 

In  the  general  case  the  coefficients  of  the  differentials  in  (6.17) 
are  evaluated  as  follows:  The  entropy  change  of  1  mole  of  an  ideal  gas 
may  be  written  in  the  following  two  forms:9 

9.  M.W.  Zemansky:  Heat  and  Thermodynamics,  McGraw-Hill  IJew  York,  1043, 
Sections  13.2  and  13.3. 
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(6. IB) 


n 


dS  =  C  -  R 

P  T  p 

dS  =  C  JT  +  J  dV.  (6.19) 

v  7  7 

Upon  elimination  of  dT  from  (6.18)  and  (6.19)  we  obtain 

# 

*  dp  =  (y  -  l)dS  -  y-  dV.  (6.20) 

P  T 

Hence  from  (6.20)  and  the  gas  law  pV  =  R7  it  f  ..lows  that 

dp  =  -  dV  +  (--f---TdS,  .6.21) 

and  therefore 

<lv>s  '  -  (6'”' 

We  wish  to  express  the  pressure  change  p'  in  terms  of  the  condensation 
s  and  the  heat  added  per  unit  volume  q.  These  quantities  are  related  to 
dV  and  dS  as  follows 


s 


dV 
V  ’ 


q 


.  TdS 
V 


Thus 


p'  =  YPqs  +  (y  -  i)q*  (6.24) 

Equation  (6.24)  must  be  used  in  place  of  (6.14)  to  obtain  the  equations 
valid  in  the  case  of  the  heated  gas  column.  Differentiating  (6.9)  we 
have 

32$  1  3p '  YP0  ?s  Y  -  1  3q 

—5-=——=  - —  t - —  .  (6.25) 

3t2  PD  at  Do  3t  PD  at 
On  combining  (6.25)  with  the  continuity  equation  (6.13),  we  get 
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i 


3  ♦  T?0  , 

—  =  —  V2*  f 

9t2  On 


1  3q 


P0  at 


(6.26) 


In  the  absence  of  neat  transfer  by  conduction  3q/3t  is  equal  to  q' ,  the 
rate  of  heat  generation  per  unit  volume.  When  thermal  conduction  is  taken 
into  account  then 


57  =  q'  ♦  KV2T, 


(6.27) 


in  accordance  with  the  heat  balance  equations  (3.1)  and  (3.2).  Thus, 
when  conduction  is  neglected,  we  have 

— ^  -  c2V2<f  =  - - —  q' ,  (6.28) 

3t2  0o 

where  c  =  /ypo/po  is  the  velocity  of  sound.  It  is  apparent  from  this 
equation  that  of  the  energy  communicated  to  the  gas  onl/  the  fraction 
Y  -  1  is  converted  into  an  acoustic  (hydrodynamic)  wave.  In  the  case  of 
diatomic  gases,  this  fraction  is  0.40.  Given  a  pulsed  laser  of  such  total 
output  that  Q  is  the  energy  incident  upon  the  test  gas  during  a  period  t, 
the  rate  of  heat  deposition  per  unit  volume  is 

q'  =  aQ/At,  (6.29) 

where  A  is  the  cross  sectional  area  of  the  laser  beam  and  a  the  absorption 
coefficient.  The  differential  equation  (6.28)  may  then  put  in  the  form 


g(x,y,z,t). 


(6.30) 


where 


E 


.  y  -  i  aQ 
p0c2  At 


f (x,y ,z,t ) . 


(6.31) 


Here  f  denotes  a  function  which  is  1  within  the  laser  beam  during  the 
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duration  of  the  pulse  and  0  at  every  other  point  in  space  or  outside  1 1 -> 
tine  interval  0  <  t  <  t . 


The  solution  of  this  differential  equation  depends  on  the  initial  i 
boundary  conditions.  We  assume  the  homogeneous  initial  conditions 
$(x,y,z,0)  =  0  and  $t(x,y,z,0)  =  0  which  correspond  to  no  d ; sturbunce  at 
t  =  0.  We  then  have  a  standard  radiation  problem  solvable  in  terms  of  t 
retarded  potentials.  In  the  absence  of  boundary  conditions;  i.e.  for  a 
homogeneous  medium  of  infinite  extent,  the  solution  of  (6.30)  for  homo¬ 
geneous  boundary  conditions  is  given  by  the  formula:1  ' 

<fr(x,y,z,t)  =  ///  gU,n,;,t-r/c)r_1  dEdndi;,  (6.3? 


wh<  re  r  =  /(x-C)2  +  (y-n)2  +  (z-O2.  The  function  g  is  a  constant  times 
a  step  function.  We  write  Eq.  (6.31)  in  the  form 

E  =  —  f (x,y  ,z,t ) ,  (<03) 

T1 


where 

Ap0c2t 

T1  (y-l)aO 


(6.34) 


is  a  constant  that  has  the  dimension  of  time.  Then 


»(x,y,z,t)  =  ^ -  HI  f(g.».n-«A»t.“il/.c.).  dEdndC.  (6.35) 


The  integration  is  to  be  extended  over  the  volume  of  space  swept  hy  the 
laser  beam.  Its  result  depends  only  on  the  beam  cross  section  and  the 
length  of  the  pulse. 


10.  A.N.  Tikhonov  and  A. A.  Samarskii,  Equations  of  Mathematical  Physics 
MacMillan,  New  York,  1963,  pp.  461-465. 


To  evaluate  the  integral  in  (6.35)  we  place  the  2-axis  along  the 
direction  of  the  propagation  of  the  laser  beam  and  we  introduce  cylin¬ 
drical  coordinates  R,  ,  and  z.  Assuming  circular  symmetry,  we  obtain 
readily 

4>(P  ,t)  =  — —  /  ■  ii -  (6.36) 

4,,T1  G  A2r(z-C,)2  ' 

where  the  region  of  integration  is  confined  to  such  values  of  c,  for 
which  0  <  t  -  r/c  <  t ,  or  equivalently 

t  -  t  <  r/c  <  t .  (6.37) 

Equation  (6.36)  is  an  apr roximat ion  valid  when  the  distance  R  from  the 
center  of  the  laser  beam  is  much  larger  than  the  diameter  of  the  laser 
beam,  i.e.  when  one  may  replace  (x-£)2  ♦  (y-n)2  by  x2  +  v2  =  R2  for 
every  point  (£,n,0  within  the  laser  beam.  We  introduce  the  new  variable 
C1  =  z-t,  then  equation  (6.36)  becomes 


$(R,t)  = 


4ht- 


‘1  G'  ^ 2  +  C 
where  the  range  of  integration  over  C ' 
(6.37).  This  consistency  means  that 

c2(t-t)2  -  R2  <  t'2  < 


.2 

is  adjusted 


c2t2  -  R2. 


(6.38) 

to  be  consistent  with 


(6.39) 


The  integration  in  (6.38)  is  symmetric  ebcut  C '  =  0,  therefore  we 


write 


VJ2 


MR.t)  =  —  ! 


dt* 


1  u:  /R2+C’5 


(0.40) 


where  and  are  non-negative  limits  consistent  with  (6.39).  We  can 
now  distinguish  three  regions: 

a)  When  R  >  ct,  the  disturbance  has  not  reached  the  observer, 
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u  =  Oand  41  =  0. 


b)  When  ct  >  R  >  c(t-T),  the  head  of  the  disturbance  has  reached 

the  observer,  but  its  tail  has  not.  In  this  case  =  0,  u  f  0 


and 


4>  - 


A  cttVc2t2-R2 

loR  - 5 - 


2ttt 


(6.41) 


c)  When  R  <  c(t-f),  the  tail  of  the  disturbance  has  begun  to  reach 
the  observer.  In  this  case  u  |  0  and  ij  }  0  and 


2itt  ■ 


{log  (ct+/c2t2-R2)  -  log  [c(t-r)  +  /c^( t-t  P-R2  J  (6.42) 


It  is  convenient  to  introduce  the  variable  R'  =  R/c  and  to  write 
4>  in  the  following  forms: 

♦<cR',t)  =  2^  l0R  ^  ^  ,  <6-43) 

for  t-T  <  R'  <  t,  and 

$(cR',t)  =  -  A-  -  {log  (t+/t2-R’2)  -  log  [t-Tt  /( t-T  )2-R ' 2]  }  (6.44) 

zTT 

for  R'  <  t  -  t.  The  graph  of  4  as  a  function  of  the  variable  t  is  sketched 
in  Figure  1  for  two  different  values  of  R'.  The  spectral  distribution  of 
the  function  $  is  clearly  not  the  same  as  that  of  the  source.  The  function 
4  is  not  a  square  wave.  Its  shape  depends  on  the  distance  R'.  It  is  ex¬ 
pected  intuitively  that  an  observer  near  the  z-axis  would  observe  something 
more  resembling  a  square  wave  than  an  observer  located  farther  away  in  a 
region  where  the  sharp  corners  of  the  excitation  function  are  washed  out. 
The  spectral  distribution  is  calculated  as  follows: 
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The  time  derivative  of  0  has  a  very  simple  form.  In  fact. 


0  for  R' 


>  t, 


_A  _L _ 

2"^  /t2-R'‘ 


for  t  >  R'  >  t-T, 


(6.45) 


2  7TT 


1  l 


A2-r'2  2 /( t- t  )2 -R ' 2 


for  t-T  >  P ' . 


The  following  equation  is  well  known: 


/  J0(bx)  sin  uxdx 
o 


0 


1 

/a  2  - 


for  a  <  b , 


for  a  >  b. 


(6.46) 


Replacing  b  by  R '  and  a  by  t  we  have 
00 

/  Jo(R'x)  sin  txdx  =  0  (6.47) 

o 

for  0  t  <  R'.  We  shall  restrict  our  observat ion:  to  the  region  P:  >  t, 
which  is  the  region  of  the  space  outside  of  the  cylinder  reached  by  the 
leading  edge  at  the  disturbance  when  the  pulse  ends.  We  integrate  (6.47) 
from  t  :  0  to  t  =  t  <  R'  and  obtain 


*J0(R'x) 

J - (1  -  cos  tx)  dx  =  0. 

o  x 


(6.48) 


Prom  (6.45)  and  (6.46)  ve  find  that  O'  may  be  written  as  follows 


OO  CO 

/  J0(R'x)  sin  txdx  -  /  JQ(R'x)  sin  (t-T)xdxl  (6.44) 

o  o  J 

Integrate  according  to  t  from  t  =  0  to  t  =  s  and  note  that  0(cR',0)  =  0. 


3$  _  A 
3t  2  TT  Ta 
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Then 


4>(cR's )  =  -A.  /“  J0(R'x)  ..  c_os_tx|£ 

2flTi  o  Lx  x  Jc 


dx 


(6.50) 


Although  J  (R'x)  £.°_s— sJi  is  not  integrable  from  0  to  <*>,  the  intepral  in 

x 

(6.50)  is  convergent  because  the  expression  in  the  brackets  is 

cos  (s-t)x  -  1  1  -  cos  tx  1  -  cos  sx 

- - -  ♦ - + - . 


and  the  product  of  each  of  these  functions  witr.  . f Q ( R  * x )  is  integrable. 
In  fact  the  middle  term  contributes  nothing  (set  eq.  6.48).  Then 


A  “  J_(  R '  x ) 

<5>(cR' ,s)  =  - -  /  -  [cos  (  s-t  )>  -  cos  sx]  dx.  (6.51) 

2n  :  x 


J  o 


Hence,  using  ar  elementary  trigonometric  identitv,  we  get 

,/  Dl  V  _  A  (  Jo(R'x)  .  TX  .  ,  T, 

$(cR' ,s)  =  - —  I  — -  sm  —  sin  (s--^x  dx. 


7T  T  2  1  X 


(6.52) 


Now  ±et  t  =  s  -  —  and  x  =  oj,  then 

x.t  di  A  (  J0(R'u)  wt  , 

<J>(cR',t+— )  1 - J  -  sin  — -  sin  hit  dui. 

2  tt  t  i  w  2 

1  o 


(6.53) 


We  have  represented  $( cR ' ,  tt^-)  as  a  Fourier  integral.  The  Fourier  sine 
transform  of  <f>  is  thus 


,(„)  =  — —  J  ( R '  uj)  Sin 
n  °  uj 


(6.54) 


Here  we  define  the  sine  traisform  by  the  relations 


F(w)  =  —  /  f(t)  sin  tit  dt,  and  f(t)  =  /  F(u>)  sin  wt  du>. 
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The  spectral  density  of  $(cR',t)  is  therefore 


♦  ( 


u.)2  =  J2(R'w)  /sin  uiT/2\2  _  A2 

°  (  T,  /  n2Ti: 


(6.55) 


It  is  shown  by  elementary  calculations  that  the  spectral  density  function 
of  a  square  wave  of  unit  amplitude  and  of  duration  t  is 


♦su,2  =(i)2(iM)2  .  <6.S6> 

Comparison  of  (6.55)  and  (6. 56)  shows  that  in  our  problem  the  spectrum  of 
the  square  wave  Rets  distorted  by  the  factor  Jo(R'w).  The  distortion  is 
small  near  the  source  where  JQ  1.  For  R 'u>  'v  2 . 4  ,  however,  JQ  =  0 .  This 
means  that  at  a  given  distance  R  the  frequency  for  which  R  =  2.4c  will  be 
missing.  Other  missing  frequencies  correspond  to  other  zeros  of  JQ.  The 
reduction  in  the  spectral  intensity  due  to  the  J0(k'w)  factor  can  be  cal¬ 
culated  from  the  following  table  which  gives  the  reduction  from  the  inten¬ 
sity  of  a  square  pulse.  The  entry  item  of  the  table  is  the  product  Rv  in 
cm  times  cycles  per  second.  Assumed  velocity  c  =  3.3  x  104  cm/sec. 


vR( cm/sec )  500  1000  1500  2000  2500  3000  3500  4000  4500 

j£  f 2irvR )  0.996  0.980  0.960  0.929  0.891  0.846  0.794  0.741  0.680 


It  is  clear  that  only  for  large  distances  or  high  frequencies  will 
the  Bessel  factor  affect  significantly  the  spectral  distribution  of  the 
energy.  The  main  determining  factor  at  small  distances  and  moderate 
frequencies  is  the  spectral  distribution  of  the  pulse.  This  distribution 
is  characterized  by  the  function 
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which  is  well  known  from  the  theory  of  diffraction.  It  is  sketched  in  F i, . 
6.2.  The  largest  part  of  the  energy  is  contained  in  the  spectrum  region 
between  0  and  wx  =  2v .  Thus,  in  terms  of  frequency,  the  interval  0  <  v  <  1/t 
contains  most  of  the  energy.  Moreover,  it  is  manifest  from  Tig.  6.2  that  the 
spectral  distribution  is  rather  flat  in  the  region  in  whi^  most  of  the 
energy  is  contained.  A  pulse  of  this  type  is  not  suitable  as  a  source  for 
a  signal  which  is  to  be  detected  by  a  narrow-band  detector. 

The  question  arises,  what  happens  when  a  pulse  of  the  type  described  is 
repeated  periodically  at  equal  time  intervals  p.  In  this  case  the  Fourier 
integral  is  replaced  by  a  Fourier  series  and  the  continuous  spectral  distri¬ 
bution  shown  in  Fig.  6,2  is  replaced  by  a  discrete  spectrum  consisting  of 
equidistant  lines  with  the  separation  Av  =  1/p.  The  relation  of  this  dis¬ 
crete  spectrum  to  the  continuous  spectrum  of  the  non-periodic  (single  pulse) 
case  is  illustrated  in  Fig.  6.3.  Again  the  energy  is  generally  broadly  dis¬ 
tributed  and  only  for  p  'v  2Tdoes  a  significant  spectral  concentration  of 
energy  take  place. 

Note:  It  is  possible  to  remove  the  restri  ion  that  the  observation 
should  be  confined  to  the  region  R’  >  x.  This  assumption  was  made  following 
equation  (6. 47).  In  the  general  case  we  cannot  assert  that  the  left  hand 
side  of  equation  (6.48)  is  zero,  but  we  can  assert  that  it  is  a  function  of 
R’  and  x  alone.  Let  it  be  denoted  by  ^(R’,x).  Then  we  can  carry  out  the 
calculations  that  lead  to  equations  (6.52)  and  (6.53)  provided  we  replace 
$(cR’,s)  by  $(cR’,s)  ♦  i|i(R’,x),  where  the  function  ip  is  independent  of  t. 
Consequently  the  spectrum  of  $(cR',t)  will  be  the  same  as  before  since  the 
time  independent  term  can  not  affect  $(w)  for  u  }  0. 
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Figure  6.1.  Distortion  of  a  square  pulse  emundtin;  :  ror.  a  line  source. 

The  function  de: ined  bv  equal'  n  (fc.‘»u)  is  s.'iown 

for  two  ui-ferent  values  or  K .  (H*  r/c.) 
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F i pure  6.3. 


Fourier  analysis  of  a  square  pulse  of  length  t 
at  intervals  p. 


repeated  periodically 
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7.  The  Acoustic  Problem  in  a  Confined  Region. 

In  the  preceeding  section  we  analyzed  the  generation  of  a  sound  pulse 
from  the  thermal  energy  deposited  in  a  gas  by  laser  pulses.  No  boundary 
conditions  were  imposed,  thus  implying  that  the  gas  was  not  confined.  If 
the  gas  were  not  confined,  or  if  the  confining  walls  were  far  removed  from 
the  path  of  the  beam,  acoustic  energy  would  radiate  in  all  directions.  The 
intensity  of  the  acoustic  signal  would  then  decrease  as  the  reciprocal  of 
the  distance  from  the  center  line  of  the  generating  beam.  The  detection  or 
measurement  of  such  a  signal  depends  to  a  large  extent  on  the  experimental 
arrangement  used  to  concentrate  a  significant  fraction  of  the  signal  gener¬ 
ated  on  a  suitable  detector  capable  of  converting  the  incident  acoustic 
signal  into  an  electric  signal.  It  would  require  a  rather  sophisticated, 
well-designed  system  to  utilize  as  much  as  20%  of  the  generated  acoustic 
energy.  Even  if  much  of  the  energy  radiated  out  is  concentrated  on  the 
detector,  by  means  of  reflectors,  for  example,  a  serious  degradation  of 
the  detector  system  will  result  from  the  fact  that  the  generated  signal 
has  a  rather  broad  spectral  distribution. 

Now  we  shall  examine  the  generation  of  sound  in  a  confined  system,  a 
cylinder  concentric  with  the  laser  beam.  It  is  expected  that  in  such  a 
system  the  efficiency  of  the  utilization  of  the  sound  generated  will  be 
higher,  and  one  might  hope  to  make  use  of  the  resonances  of  the  confining 
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tube  to  enhance  detection.  The  following  questions  are  worth  examining: 

1.  How  does  the  confining  tube  affect  the  generation  and  propagation 
of  the  sound  wave? 

2.  Is  it  possible  to  enhance  the  measurement  of  the  sound  wave  gener¬ 
ated  by  proper  choice  of  the  geometry  of  the  confining  vessel? 

The  starting  point  of  the  analysis  is  the  differential  equation  (6.30) 

V2$  -  k  lli  =  -  g(x,y,z,t)  (7.1) 

c2  at2 

with  the  initial  conditions 

4>(x,y,z,0)  =  0;  $t(x,y,z,G)  =  0.  (7.2) 

The  meaning  of  the  variables  was  explained  in  section  ft,  where  it  was  also 
noted  that  g  is  zero  except  in  the  part  of  the  space  swept  out  by  the  laser 
beam.  In  this  cylindrical  region,  g  has  a  constant  value  during  the  time 
interval  0  <  t  <  t;  it  is  0  when  t  >  t.  We  shall  simplify  the  computations 
by  normalizing  g  so  that  is  has  the  values  0  and  1.  Since  g  has  the  dimen¬ 
sion  of  reciprocal  time,  it  is  more  correct  to  say  that  we  set  the  peak 
value  of  g  at  1  sec-1. 

The  presence  of  the  walls  of  the  tube  requires  the  introduction  of  a 
boundary  condition  which  expresses  the  fact  the  normal  component  of  the 
fluid  velocity  at  the  walls  is  zero.  This  means 

n-V$  =  0.  (7.3) 

at  all  boundaries.  We  shall  seek  a  solution  of  (7.1)  in  the  region 
x2  +  y2  $  b2  for  t  >  0.  A  useful  technique  for  the  solution  of  such  a 
problem  is  the  eigenfunction  expansion.  This  technique  relies  on  the  deter¬ 
mination  first  of  the  simple  harmonic  solutions  of  the  homogeneous  equation 
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obtained  from  (7.1)  with  g  =  0.  It  is  well  known  that,  if  we  set 

$  =  e-iut  ij/(x,y,z),  (7.4) 

then  fp  satisfies  the  homogeneous  Helmholtz  equation 

V2ip  +  k'<|»  =  0  (7.5) 

which  is  obtained  from  (7.1)  by  setting  g  =  0,  substituting  (7.4)  and 
k  =  u/c.  It  is  also  known  that,  given  a  finite  region  Si  and  given 
homogeneous  boundary  conditions,  non-vanishing  solutions  of  (7.4)  are 
possible  only  for  a  discrete  set  of  values  of  k,  called  the  eigenvalues 
of  (7.4)  and  that  solutions  corresponding  to  different  eigenvalues  are 
orthogonal  to  each  other  in  the  sense 

/  i p  <p  *  dv  =  0,  (7.6) 

‘  n  m 

whenever  kn  p  km.  Following  standard  practices  in  mathematical  physics, 
we  obtain  an  ascending,  infinite  sequence  of  eigenvalues 

k  <  k_  <  k,  <  • • • »  (7.7) 

X  —  l  —  0  — 

and  a  corresponding  sequence  of  orthogonal  and  dimensionless  eigenfunctions 

so  normalized  that  they  satisfy  the  relations 

/  4^  ip  j  dv  =  6^.  V,  (7.8) 

where  dv  =  dx  dy  dz  and  integration  is  over  the  entire  volume  V  of  region 
.  It  is  known  that  the  sequence  of  such  functions  is  complete  for  equa¬ 
tions  of  the  type  we  are  considering  and  that  a  well-behaved  function  of 
x,  y,  z  is  expandable  in  the  region  51  in  terms  of  these  eigenfunctions. 
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Convergence  in  the  mean  is  assured  by  the  standard  theorems  of  functional 
analysis . 


With  these  preliminaries  we  can  reduce  the  solution  of  our  original 
problem  to  the  calculation  of  the  excitation  of  an  infinite  number  of 


oscillators  by  a  given  exciting  force.  We  write 

<J>  =  l  un(t)«|>n(x,y,z),  (7.9) 

g  =  l  Rn(t)^n(x,y,z).  (7.10) 

The  coefficients  u  (t)  and  g  (t)  are  calculated  as  follows: 
n  ^n 

v  un  =  /  Wn*  dv,  (7.11) 

V  gn  =  /  R<J>n*  dv-  (7.12) 

The  differential  equation  (7.1)  then  requires  that 

V  ♦  un2  yn  =  Rn  (7.13) 


for  n  =  1,2,3,  **’.  Since  the  r,n's  are  calculable  from  (7.12),  the  solu¬ 
tion  of  the  complete  problem  now  reduces  to  the  determination  of  the  ex¬ 
citations  of  the  harmonic  oscillators  whose  equations  are  given  above. 

The  excitations  depend  on  the  characteristic  frequencies  of  these  oscillators, 
<i)n,  and  on  the  functions  gn(t). 

The  total  energy  fed  into  the  n-th  oscillator  is  proportional  to 

/  |gn(t)|2dt 
o 

and  to  another  factor  which  represents  the  response  of  the  n-th  oscillator 
to  a  normalized  exciting  force.  This  response  factor  depends  on  u>n  and  on 
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the  time  dependence  of  the  exciting  force  gn.  (It  is  to  be  noted  that  gn 

has  the  sane  dimension  as  g,  which  is  not  that  of  force.  At  the  same  time, 
the  dimension  of  un  is  that  of  $.) 

In  principle  the  calculations  require  the  solution  of  infinitely  many 
differential  equations  of  the  form  (7.13)  for  very  general  functions  gn(t) 

on  their  right  side.  Practically  the  situation  is  much  simplified  because 
the  physically  relevant  part  of  the  solution  may  be  obtained  by  solving 
only  one  or  a  few  equations  of  the  type  described.  There  are  two  main 
reasons  for  this  simplification.  First,  in  a  situation  with  some  degree 
of  symmetry  many  of  the  functions  gn(t)  are  identically  zero,  because  the 
spatial  part  of  the  excitation  function  is  orthogonal  to  many  ^m's.  Second, 
the  sequence  kj_,  k2,  increases  rather  rapidly  and  consequently  only  a 
few  of  the  harmonic  oscillators  which  are  described  by  equation  (7.13)  have 
resonant  frequencies  in  the  region  in  which  f  has  a  significant  harmonic 
component . 

We  shall  return  to  the  theory  after  exhibiting  some  of  the  well-known 
eigenvalues  and  eigenfunctions  appropriate  for  sound  propagation  in  a  long 
circular  cylinder. 

The  eigenfunctions  appropriate  for  a  right  circular  cylinder  of 
length  L  and  radius  b  are  the  non-vanishing  solutions  of  equation  (7.5) 
which  satisfy  the  boundary  conditions 

P  =  0  for  r  =  b,  and  P  =  0  for  z  =  +  L/2. 

3r  3z  - 

After  separating  equation  (7.5)  in  cylindrical  coordinates,  one  obtains 
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solutions  of  the  form 

4/  =  J  ( A2-h2r)  cos  hz  (7.14) 

s  n 

and 

ip  =  J  ( A2-’n2r)  sin  hz,  (7.15) 

3  n 

where  the  subscript  s  denotes  a  function  symmetric  for  the  reversal 
of  the  z-axis,  while  subscript  a  denotes  an  antisymmetric  function.  Jn  is 
the  Bessel  function  of  order  n,  where  n  must  be  an  integer  to  insure  unique¬ 
ness  of  the  solution.  These  functions  are  not  yet  normalized.  We  note  that 
if  the  cylinder  contained  one  or  more  septum;  i.e.  a  boundary  plane  contain¬ 
ing  the  cylinder  axis,  then  n  need  not  be  an  integer.  We  shall  avoid  this 
situation. 

The  possible  values  of  h  are  determined  by  the  boundary  conditions 
at  the  cylinder  ends.  These  require  that,  for  symmetric  functions  sin 
hL/2  =  0;  i.e. 

h£  =  2£n/L;  l  =  0,1,2,-.. .  (7.16) 

For  antisymmetric  functions  cos  hL/2  =  0;  i.e. 

h£  =  (21  *  1 )  7T /L ;  l  =  0,1. 2,-...  (7.17) 

Boundary  conditions  on  the  cylinder  surface  r  =  b  require  that 

Jn'  (  A^-h2b)  =  0.  (7.18) 

Therefore  k  can  only  have  such  values  for  which  the  number  x  =  A2-h2b 
is  a  solution  of  the  equation  Jn'(x)  =  0.  For  each  index  n,  the  deriva¬ 
tive  of  the  n-th  order  Bessel  function  has  an  infinite  number  of  discrete 
positive  roots.  We  write  these  in  ascending  order  as  x^;  m  =  1,2,,,‘. 

The  eigenvalues  of  the  differential  equation  (7.5)  consistent  with  the 
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boundary  conditions  described  above  are  calculated  from  the  equations 

(knml  •  h!)b2  1  *nm;  n  =  "  *  U2 ,  ■  ■  ■ ,  (7. ID 


where  h^  are  the  numbers  determined  from  (7.16)  and  (7.17).  The  symmetric 
eigenfunctions  are 

.in8 


*snmi  =  e  Jn(*nmr/b)  cos  h»z’ 


(7.20) 


while  the  antisymmetric  eigenfunctions  are 


4*  .  =  e 

anmi 


in8 


J  <X  r/b)  sin  h.z. 
n  nm  i 


(7.21) 


Of  the  triple  infinity  of  eigenfunctions  (or  oscillation  modes)  given  in 
(7.20)  and  (7.21),  certain  types  are  of  particular  interest  to  us.  Only 
l  =  0  produces  z-independent  eigenfunctions.  These  are  called  transverse 
modes.  Their  equations  are 

4-  =  ein0  J  (x  r/b).  (7.22) 

snm  n  "nm 

Modes  which  have  rotational  symmetry  about  the  z-axis  are  all  of  the  form 


and 


*som«.  =  Jo(*omr/b)  cos 


ip  .  =  J  (y  r/u)  sin  - 
aomi  o  Aom  2L 


(7.23) 


(7.24) 


The  first  few  roo”s  of  JQ'(r)  =  0  are  3.832,  7.016,  10.173,  13.324. 


Those  of  J^'(r)  =  C  are  1.841,  5.330,  etc 
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Let  us  now  calculate  the  excitation  of  a  cylindrical  resonator  with 
a  force  function  that  depends  only  on  r  and  t;  i.e.  which  is  independent 


11.  See  Ref .  8 
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of  0  and  z.  In  this  case  every  integral  of  the  type 
b  2  it 

/  /  g(r,t)  JR(  x^r/bJe^H®  rdrde 

vanishes  for  n  ^  0,  leaving  only  the  eigenfunctions  with  n  =  0  to  interact 
with  g.  Similarly  the  integrals  over  z  involving  the  product  of  g  with 
sin  h^z  and  cos  h^z  will  vanish  except  for  cos  Oz.  Thus,  when  evaluating 
the  expansion  of  g(r,t)  in  terms  of  the  eigenfunctions  of  the  problem, 
only  the  symmetric  eigenfunctions  with  n  =  0  and  l  -  0  need  to  be  consid¬ 
ered.  The  relevant  eigenfunctions  are  JQ(x0rnr/b).  Their  normalization 
factor  is  calculated  from  the  standard  formula12 

/  lJo(W'/b)]2rdr  2  7  b2[Jo(*o»,)]2-  (7-J5) 


valid  for  every  xQm  f°r  which  Jo'^Xom)  =  0.  The  normalized  functions  are 

*„  -  Jo<W'b)/J.<*OT,>.  (7-J6> 

When  the  excitation  is  of  unit  magnitude  in  a  small  cylinder  of  radius  a 
centered  around  the  polar  axis  and  lasts  from  t  =  0  to  t  =  t,  we  have  from 
(7.12)  and  (7.26) 

Vgm(t)  =  j'  Ty~L" ')'  /  J  (xomr/b)rdrf(t),  (7.27) 

^o'^om'  o 

where  f(l)  =  1  sec-1  for  0  <  t  <  t,  and  0  for  t  >  t.  For  small  values  of 
a  we  may  replace  Jo^X0mr/b)  by  1  and  obtain  the  approximation 


„  (tl  .  v  f(t) 
'  v 


(7.28) 


where  v  =  a27»L  is  the  volume  swept  out  by  the  laser  beam. 

( 


12.  Ref.  8,  p.91. 
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The  next  problem  is  to  deal  with  a  set  of  harmonic  oscillators  with 
pulse  excitation.  Before  losses  in  the  oscillators  are  considered,  their 
equations  ha’'e  the  form 

where  gm  is  given  in  (7.28).  The  characteristic  frequencies  are 

um  =  Xomc/b.  (7.30) 

Taking  c  =  3.3  x  104  cm/sec,  we  get  for  the  lowest  symmetric  transverse 
frequency  (corresponding  to  =  3.832) 

.  2.0xl04 

V1 - E - 

and  for  the  next  lowest 

„  _  3.66xl04  , 

2  "  b 

where  b  is  in  cm,  v  in  sec"J.  Thus  the  lowest  such  frequency  for  a  tube 
10  cm  in  diameter  is  4000  sec"1,  the  next  7330  sec"1.  A  small  tube  of  1 
cm  diameter  would  require  acoustic  frequencies  which  are  prohibitively 
high. 


Generally  the  oscillators  will  have  a  certain  loss  rate  associated 
with  them.  Losses  are  not  only  unavoidable  but  necessary  when  the  res¬ 
onator  is  a  part  of  a  system  used  for  measurement,  because  measurement 
consists  of  converting  a  part  of  the  acoustic  disturbance  into  an  elec¬ 
tric  signal.  Thus  a  resonator  will  have  an  output  which  constitutes  the 
measured  signal,  or  else  the  energy  is  communicated  to  another  frequency- 
sensitive  device,  the  detector,  which  converts  it  to  an  electric  signal. 
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In  any  case,  we  are  interested  in  devices  which  subject  the  disturbance 
to  a  harmonic  analysis.  The  noise  discriminating  property  of  the  detecting 
system  improves  as  its  bandwidtn  is  restricted.  Unfortunately,  as  we  shall 
see,  narrowing  the  bandwidth  drastically  reduces  the  signal. 

for  the  sake  of  simplicity,  let  us  assume  that  the  detector  is  the 
cavity  itself,  that  is,  we  do  away  mentally  with  a  second  tuned  device  to 
which  the  cavity  may  be  coupled.  If  there  is  such  a  separate  device,  the 
detected  signal  will  be  less  because  some  of  the  energy  will  not  be  trans¬ 
ferred  and  new  sources  of  noise  will  be  introduced.  We  thus  have  a  sequence 
of  oscillators  with  an  output  and  we  concentrate  on  the  one  with  the  lowest 
frequency.  When  the  effects  of  oscillator  loss  (energy  conversion)  are  in¬ 
corporated  in  the  equations  of  the  harmonic  oscillators,  we  obtain  for  the 
oscillator  of  the  lowest  frequency13 

y"  +  Au'  +  oj^u  =  g}(t).  (7.31) 

The  quantity  A  has  the  dinension  of  reciprocal  time.  Let  W  denote  the 
energy  stored  in  the  free  running  harmonic  oscillator.  It  is  known  that 

A  =  ~  W- 1  (7.32) 

at 

/ 

One  can  readily  show  the  following:  When  g  is  a  sinusoidal  function 
of  time,  the  excitation  of  the  oscillator  is  significant  only  if  the  fre¬ 
quency  of  the  exciting  force  is  nearly  equal  to  that  of  the  free-running 
oscillator.  Then  by  means  of  Fourier  analysis  one  concludes  that  in  the 
case  of  a  general  time -dependent  force,  the  excitation  of  the  oscillator 

lX  P.M.  Morse,  Vibration  and  Sound,  McGraw-Hill,  New  York,  1948,  Ch .  II. 
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depends  largely  on  the  harmonic  component  of  the  exciting  force  which  lies 
near  v^.  After  a  short  calculation  one  finds  that  the  average  energy  of 
the  oscillator  is  approximately 


W  % 


(7.33) 


where  G(v)  is  the  Fourier  transform  of  the  exciting  function  g(t)  and 
V1  =  The  power  output  of  such  an  oscillator  is 

^  =  AW  i-|G(  vx)  |2  (7.34) 


Only  part  of  this  is  useful  power  output  in  the  detector,  hut  this  brings 
in  only  a  factor  of  about  one-half,  which  is  not  significant,  "he  important 
matter  is  that  the  rate  at  which  power  is  absorbed  by  the  oscillators  is 
proportional  to  the  Fourier  transform  of  the  excitation  function  at  the 
resonant  frequencies  of  the  oscillators. 


The  effect  of  the  confining  tube  is  that  the  acoustic  signal  generated 
within  it  is  resolved  into  its  spectral  components  according  to  the  resonant 
frequencies  of  the  cavity.  As  we  have  seen  in  section  6,  the  power  spec¬ 
trum  of  a  single  pulse  of  duration  t  is  rather  broad,  most  of  the  energy  is 
in  the  spectral  interval  0  to  1/x.  Periodic  repetition  of  the  pulse  pro¬ 
duces  a  spectral  distribution  of  the  acoustic  energy  among  the  harmonics  of 
the  repetition  frequency  in  the  manner  illustrated  in  Fig.  6.3. 


In  order  to  obtain  good  utilization  of  the  acoustic  energy  generated, 
the  experiment  should  be  so  designed  th3t  a  large  fraction  of  the  acoustic 
energy  generated  by  the  source  is  concentrated  in  the  frequency  region  cap¬ 
able  of  exciting  a  single  selected  cavity  resonance  because  the  detector 
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will  only  be  coupled  to  one  resonance.  Such  concentration  of  the  acoustic 
energy  is  not  easily  achieved.  With  a  single  pulse  it  is  just  about  hope¬ 
less  to  obtain  success  in  this  manner.  With  pulses  repeated  at  the  rate  of 
a  resonant  frequency  of  the  cavity,  the  situation  is  better,  but  reason¬ 
able  cavity  dimensions  require  rather  high  pulse  repetition  frequencies. 

We  calculated  that  the  lowest  resonance  of  a  cylindrical  cavity  of  10  cm 
diameter  is  at  4000  Hz,  therefore  we  would  have  to  operate  with  the  laser 
pulses  shorter  than  0.1  msec  repeated  4000  times  per  second.  The  success 
of  the  scheme  of  this  type  depends  not  only  on  the  construction  of  a  power¬ 
ful  laser  with  appropriate  pulsing  mechanism  but  also  on  the  construction 
of  a  nighly  sensitive  acoustic  detector  coupled  with  the  proper  mode  of 
oscillation  of  a  cylinder. 

8.  Interferometer  Method. 

It  is  plausible  to  contemplate  a  two  laser  experiment  in  which  one 
high  intensity  laser  pulse  heats  the  air  and  the  second  probe  laser  beam 
travels  the  same  path  which  is  a  branch  of  an  interferometer. 14  Then  the 
heating  effect  produced  by  the  intense  laser  pulse  would  cause  a  shift  in 
the  interference  pattern  of  the  probe  beam  due  to  a  change  in  the  optical 
path  length  of  one  of  the  branches  of  the  interferometer.  See  Figure  8.1. 
The  absorption  chamber  nay  be  evacuated  to  determine  the  effect  of  the  in¬ 
tense  laser  pulse  on  the  end  windows;  or,  if  measurements  are  to  be  made 
in  the  atmosphere,  the  end  windows  may  be  removed. 

Ii  Q  is  the  energy  transmitted  in  the  intense  pulse  of  cross  sectional 

TIT  Longaker  and  Litvak:  Refractive  Index  Changes  in  Absorbing  Media  by 
a  Pulsed  Laser  Beam.  Bull.  Am.  Phys .  Soc .  11,  No.  1,  129,  1966. 
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area  A,  then  the  energy  balance  over  a  unit  length  of  the  gas  traversed 
by  the  intense  pulse  gives: 


Qo  s  ApJ'cpdT, 

where  c  is  the  specific  heat  of  the  gas  at  constant  pressure 

r 

J'  =  4.19  joules/cal. 

6T  is  the  increase  in  temperature  of  the  gas. 

Hence 


0 


A  J  'pc. 


The  temperature  increase  6T  produces  n  density  change  dp.  In  the  case  of 
an  ideal  gas,  where  the  temperature  change  is  assumed  to  take  place  quickly 
enough  to  occur  at  constant  pressure. 


T  +  p 


0. 


The  change  of  density  causes  a  change  in  the  refractive  index  according  to 
the  relation. 


n-1  _  p 

"o'1  ’  Po  ’ 

where  n  is  the  refractive  index  of  the  gas  at  the  frequency  of  the  probe 
laser  beam. 

From  these  relations  it  follows  that 

dn  »  (n-1)  &  =  -  (n-1)  il  , 

P  * 


or 


5n  =  - 


(n-1) 

T 


0°  . 

AJ'pcp 


55  - 


Now  if  L  is  the  length  of  the  absorption  chamber,  the  change  in  optical 
path  length  for  the  probe,  or  monitor,  laser  beam  is  L  fin,  which  will  be  a 
certain  fraction,  F,  of  its  wave  length.  Hence, 

L  fin  =  FX, 


and  one  has 


F 


-  X 


Q 

AJ'pepT  ‘ 


The  minimum  value  of  a  which  can  be  detected  depends  upon  how  small  a  value 
of  F  one  can  measure  for  a  given  value  of  the  other  physical  quantities  in 
the  experiment. 


For  example,  assuming  Q/A  =  20  joules/cm2,  L  =  100  cm,  X  =  5  x  10"5  cm, 

(n-1)  =  3  x  10-4,  T  =  300°K,  p  *  1.3  x  10-3  gm/cm2 ,  cp  =  0.24  cal/gm  °C, 

one  has  F  *  3  x  10*4  a. 

If  one  can  measure  a  path  length  change,  or  phase  shift,  of  a  thousandth  of 
a  wave  length  (F  *  0.001),  one  can  then  detect  a  of  3  x  10"®  cm-1.  One  can 
increase  the  energy  and  path  length  to  perhaps  detect  an  a  of  an  order  of 

magnitude  less.  But  a  phase  change  of  a  thousandth  of  a  wave  length  is  dif¬ 

ficult  to  measure  by  methods  usually  associated  with  interferometry. 


In  order  to  measure  extremely  small  phase  shifts,  of  the  order  of 
X/1000,  or  possibly  smaller,  the  second  branch  of  the  interferometer  probe 
laser  beam  is  also  varied  in  optical  path  length,  but  at  a  frequency  of  the 
order  of  105  Hz.15  This  may  be  accomplished  by  the  insertion  of  a  crystal¬ 
line  Kerr  Cell  in  the  optical  path,  or  by  making  one  mirror  a  Piezo-electric 
crystal,  either  one  driven  by  a  105  Hz.  oscillator.  The  interference  pattern 

1ST  P.H.  Lee  and  Skolnick:  "A  Phase  Comparison  Optical  Discriminator", 

IEEE  J.  Quantum  Electronics,  QE-2,  No.  12,  p.  784-785,  1966. 


-  56  - 


is  detected  by  a  photomultiplier  tube  whose  output  will  be  the  10s  Hz. 
signal  caused  by  the  driven  crystal.  This  is  fed  into  a  synchronous  de¬ 
tector,  together  with  a  signal  directly  from  the  105  Hz.  oscillator. 

Now  a  servo-mechanism  is  used  to  apply  a  variaoie  high  D-C  voltage  to  the 
crystal  in  order  to  further  vary  the  optical  path  of  the  second  branch  in 
the  saute  amount  that  the  intense  laser  pulse  varies  the  optical  path  of 
the  first  branch,  thus  keeping  the  time  average  optical  path  of  the  second 
branch  equal  to  the  optical  path  of  the  first.  See  Figure  8.2.  The  inter¬ 
ference  pattern  is  then  kept  in  phase,  ar.d  the  D-C  servo-voltage  required 
to  do  this  is  proportional  to  the  phase  shift  that  would  otherwise  occur. 
Hence  the  small  phase  shift  is  converted  into  a  measurable  voltage. 

The  limitations  here  are  determined  by  the  constancy  in  optical  path 
one  can  achieve  by  carefully  constructing  the  interferometer  and  by  the 
homogeneity  of  the  wave  fronts  in  the  first  branch  through  which  the  in¬ 
tense  laser  is  pulsed.  A  detectivity  of  o  down  to  the  order  of  10“9/cm 
should  be  achievable  under  favorable  conditions.  This  appears  as  the  most 
promising  alternate  to  the  spectrophone,  especially  for  intense  laser 
pulses  at  frequencies  where  scattering  is  appreciable.  Smaller  values  of 
a  might  be  detectable  by  measuring  effects  over  long  ranges,  as  might  be 
possible  in  a  mile  long  tube. 
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Probe  Laser 


// 


Figure  8.2.  Interferometer  With  Phase  Comparison  Optics15 


9.  The  Absorptivity  Spectrophone  and  Measurement  of  a. 

The  absorptivity  spectrophone,  a  practical,  direct  method  for  precise 
measurement  of  weak  optical  absorptivity  in  gases,  especially  for  wave 
lengths  at  which  scattering  is  not  significant,  has  been  designed  by  Edwin 
L.  Kerr  and  John  G.  Atwood.16  They  used  the  spectrophone  to  measure  the 
absorption  coefficient  of  the  9.6  micron  line  of  a  CW  CO2  laser  in  mixtures 
of  C02  -  Nj.  At  300  ppm  C02  in  Nj,  the  standard  concentration  of  C02  in 
air  at  sea  level,  the  absorptivity  reported  was  5.8  x  10"7  cm"1.^16^ 

The  same  equipment,  the  theory  of  which  has  been  discussed  in  Section  3, 
was  used  by  us  in  an  attempt  to  measure  the  absorption  coefficient  of  the 
10.6  micron  line  of  a  CW  CO2  laser  in  the  actual  atmosphere.  At  24°  +  1°  C 
with  a  relative  humidity  of  (70  +  5)%  and  an  unknown  amount  of  air  pollu¬ 
tants  present  in  the  air  in  the  Western  San  Fernando  Valley,  the  measured 
value  of  the  atmospheric  absorption  coefficient  is  (1.0  +0.2)  x  10"6  cm-1. 
Time  and  experLnental  difficulties  did  not  permit  the  accomplishment  of  a 
more  ambitious  goal,  that  of  measuring  a  in  controlled,  known  mixtures  of 
n2  "  °2  "  co2  "  h2°  while  simultaneously  monitoring  the  spectral  distribu¬ 
tion  of  the  laser  beam. 

The  spectrophone  measures  absorptivity  by  sensing  the  thermal  expan¬ 
sion  of  a  confined  sample  of  gas.  A  complete  description  of  the  device  is 
contained  in  the  Perkin-Elmer  Report^ 16^.  Only  a  brief  description  is  con¬ 
tained  here.  The  laser  was  made  in  house,  except  for  the  pyrex  tube  which 
was  loaned  to  us  by  a  prominent  manufacturer  of  same.  The  length  is  0.7 
meters  with  an  inside  diameter  of  0.5  inches.  Using  one  plane  and  one 

HT!  Earl  L.  Sloan  and  Edwin  L.  Kerr,  Perkin-Elmer  Corporation  Report  No. 
8884,  Absorption  of  Light,  31  July  1967. 
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spherical  mirror,  it  supplied  an  output  of  just  under  one  watt  during  the 
experiments  while  excited  by  a  D>C  voltage  of  12,000  and  containing  the 
following  (closed)  gas  mixture: 

1.0  Ton’  C02,  2.5  Torr  N2,  11.0  Torr  He,  0.2  Torr  H20. 

The  spectral  distribution  is  dependent  upon  many  factors,  including  proper 
alignment  of  the  end  mirrors.  Spectral  measurements  were  not  made  simul¬ 
taneously  with  atmospheric  absorption.  Hence  all  that  can  be  said  is  that 
measurements  made  at  different  times,  when  the  laser  is  tuned  for  maximum 
power,  as  during  the  experiments,  indicate  that  the  laser  power  is  almost 
entirely  at  10.6  microns.  Laser  power  was  measured  with  a  Korad  Model  K-PM 
Power  Meter. 

The  absorptivity  measurements  were  made  by  detecting  the  pressure  ?ise 
in  a  gas  sample  chamber  while  the  CW  C02  laser  beam  is  passed  through.  The 
sample  chamber  and  pressure  transducer  were  obtained  from  the  Pcrkin-Elmer 
group,  where  they  were  used  for  similar  measurements.  The  sample  chamber 
is  aluminum,  2  cm  I.D.,  10  cm  O.D.,  20  cm  path  length,  4.2  kg  mass,  and 
has  6mm  KBr  windows  at  Brewster's  angle.  The  chamber  is  massive  to  insure 
a  uniform  constant  temperature.  Thermal  and  acoustic  shielding  is  provided 
by  an  outer  cylindrical  aluminum  tank  which  is  sealed.  The  laser  beam  power 
was  monitored  by  a  Korad  K-PM  power  meter.  The  pressure  rise  was  measured 
by  an  MKS  Baratron  pressure  transducer  with  a  minimum  resolvable  differential 
of  10~5  mm  Hg.  or  10~8  atm.  A  steady  state  temperature  distribution  is 
reached  in  a  few  tenths  of  a  second  after  the  shutter  open-}  to  allow  the  beam 
to  shine  through  the  chamber. 

As  derived  in  Section  3,  with  the  simplifying  assumption  of  a  perfectly 
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defined  cylindrical  beam,  the  pressure  rise  is: 


dP 


Pq*JW  /  _a£ 
T04*K  (  2b2 


(3.10) 


where  PQ  *  initial  pressure 

Tq  *  initial  temperature 
a  =  linear  absorption  coefficient 
W  *  beam  power 

K  a  thermal  conductivity  of  gas  in  sample  chamber 
J  a  0.239  cal/ joule 
a  =  beam  radius 
b  a  sample  chamber  radius . 

In  actual  practice,  the  power  density  of  the  laser  beam  is  approximately  a 
Gaussian  distribution,  with  the  intensity  given  b”: 


I(P)  = 


W  e-p2/w2 


4irw2 


where  w  =  half  beam  width  at  e"1  points,  p  =  distance  from  axis. 

2 

refinement,  the  quantity  in  parenthesis  (1  -  _iL)  is  replaced  by 

2b2 


With  this 
16 


(1 . 4  ♦  (4  -  T-)  .-b2/“2]  *  (1  -  4>  for  »  «  b. 

b7  b2  2W7  b2 

The  reproducibility  of  measurement  of  the  pressure  rise,  dP,  was  of  the 

order  of  10%,  hence  a  precise  measurement  of  w  is  unnecessary,  as 
2 

(1  -  iL.)  a  0.98  +  0.01  in  our  experiment. 
b2 

The  absorption  chamber  was  baked  out  at  temperatures  in  excess  of  100°C 
purged  with  dry  nitrogen  several  times  to  attempt  to  rid  the  KBr  windows  of 
contaminants,  but  enough  remained  to  give  an  equivalent  absorptivity  of 
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10- 7  cm*1  with  dry  Nj  at  atmospheric  pressure  in  the  cell.  The  was 
evacuated  and  air  was  admitted  to  the  cell,  and  absorptivity  measurements 
were  made  on  several  days  —  all  of  which  were  midly  smoggy.  With  the 
consistent  weather  in  the  Western  San  Fernando  Valley,  the  climate  con¬ 
ditions  were  similar  on  all  days.  Temperature  (24  ♦  l)°C.  Pressure 
740  +_  5  mm  Hg.  Relative  humidity  (70  5)%.  Very  mild  smog  conditions. 

Under  these  conditions  the  measured  value  of  a  =  (1.0  ♦  0.2)  x  10*6  cm*1. 

10.  Conclusions. 

The  following  conclusions  were  reached  concerning  the  measurement  of 
the  absorption  of  raaiation  in  air  for  wave  length  regions  where  a  <  10~8 
cm”1 : 

1.  It  is  probably  possible  to  make  absorption  measurements  in  air 
using  a  steady  laser  beam  and  noting  the  static  pressure  in¬ 
crease  in  a  sealed  tube  if  the  following  conditions  can  be  met: 

a)  The  tube  i3  long  enough  to  render  the  end  effects  negligible. 

b)  The  heating  effect  of  the  scattered  light  is  eliminated. 

c)  The  temperature  of  the  confining  tube  is  maintained  constant 
to  within  lO’4  °C. 

2.  The  static  pressure  rise  due  to  a  single  laser  pulse  can  probably 
not  be  measured  if  a  is  as  low  as  10" 9  cm*1. 

3.  The  measurement  of  the  acoustic  wave  generated  by  a  single  laser 
pulse  depends  on  th<  issing  of  the  acoustic  wave  and  on  an 
extremely  sensitive  broad-band  detector.  It  is  doubtful  that 
this  can  be  achieved. 


-  63 


4.  The  measurement  of  the  acoustic  signal  generated  by  periodically 
repeated  pulses,  whose  repetition  frequency  is  matched  to  that 

of  •  resonant  cavity,  is  a  possibility  tha'  should  be  investigated 

f urther . 

5.  By  use  of  a  phase  comparison  optical  discriminator  described  in 
Section  8,  values  of  a  in  the  range  10-8  to  10-9  cm  probably  can 
be  measured  by  interferometric  methods. 

The  situation  in  the  wave  length  region  around  10  microns  is  quite  different 
from  that  investigated  here.  First,  Rayleigh  scattering  is  no  longer  a 
limiting  factor.  Second,  the  absorption  in  air  is  considerably  higher  at 
10  microns  than  it  is  in  certain  parts  of  the  visibl-  region.  Hence  several 
of  the  contemplated  methods  would  work.  The  absorptivity  spectrophone  was 
used  to  measure  a  in  air  for  the  10. 6u  C02  laser  line,  a  =  (1.0  +  0.2) 
x  10-6  cm-1. 
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